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1. Introduction 

In this introductory section we want to recall some well-known estimates for the 
wave equation and their relation to Fourier integrals and restriction theorems. This 
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shall serve as motivation to study similar estimates for perturbations of the wave 
equation by lower order terms and outline strategies to be pursued. 

1.1. Energy and dispersive estimates 

We will recall some well-known properties of solutions to the Cauchy problem for 
the wave equation 

utt-Au = 0, u{0,-)^uo, ut(0, •) = ui (1.1) 

in M X E" . We assume for simplicity that data belong to C^ (M" ) . A simple integra- 
tion by parts argument allows to show that solutions possess a cone of dependence 
property and are therefore also compactly supported for all fixed times t. We will 
not make this argument rigorous, but point out that by the same reason the total 
energy defined as 

E{u;t)^^J{\Vu\^ + \ut\')dx (1.2) 

is preserved, 

— E(m; t) =Re [utUu + Vu • Vut)dx = Re / ut{utt - AM)dx = 0. (1.3) 

This preserved energy is spread out (uniformly) over a region increasing in 
time. Dispersive estimates due to von Wahl [SQj , Brenner [S] , Pecher [37] and more 
generally Strichartz estimates [53] describe this effect. See also Keel-Tao [13] for 
the precise interrelation between Strichartz and dispersive estimates. 

We denote by LP(M") the Bessel potential space of order r over LP(M") 
defined in terms of the Fourier transform as 

LP(M") = (D^)-''iP(M") = .^-^[{S)-''.^LP{W)]. (1.4) 

Then the following statement contains the essence of the above cited papers. Here 
and later, we denote (0 = (1 + \i?Y''^- 



Theorem 1.1. Solutions to (1.1) satisfy the dispersive estimates 

\\u{t, OIIl- < Crt-'-^ {\\uq\\li + II"i||li_,) , (1-5) 

||Vu(i, ■)\\l^ + \\ut{t, Oik- < Crt-"^ (|I"o||li^, + IImiIIli) , (1-6) 

forr> ^. 

Dispersive estimates are much harder to prove than energy estimates as they 
encode structural information about the representation of solutions. We will collect 
some of the crucial ingredients of their proof in this introductory section. 
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1.2. Equations with constant coefficients 

Of particular importance for us are equations with constant coefficients as they 
can be solved explicitly in terms of the Fourier transform and serve as important 
model examples. We keep it as simply and sketchy as possible to not extend the 
exposition too much. Later we will discuss variable coefhcient versions of the wave 
equation 

uu - A« = 0, (1.7) 

the damped wave equation (or dissipative wave equation) 

wtt-AM + ut = 0, (1.8) 

and the Klein-Gordon equation 

Mft-Au + u = 0. (1.9) 

Natural questions to ask for these problems are: 

1. What are reasonable energies to estimate? While for the damped wave equa- 
tion E(u;i) is a good candidate, we need to include u itself into the energy 
term for the Klein-Gordon equation. Then the Klein-Gordon energy is pre- 
served, while for the damped wave equation we obtain dtR(u; i) < and the 
(wave) energy decreases. 

2. What are sharp (decay) rates for the energy? This is of interest for the 
damped wave equation and the answer is by no means trivial. It heavily 
depends on the assumptions made for data, see e.g. the work of Ikehata [T7] . 
[18] . [20] for some examples of this. The obvious estimate 

E(u;i) <E(m;0) (1.10) 

is sharp for the damped wave equation with arbitrary data from the energy 
space. On the other hand there is the estimate of Matsumura [28] 

E{u-t) < (1 + f)-i-t (ho||?finLi + hilli^niO ■ (l-H) 

3. How do dispersive estimates for the above models look like? Here we will 
just point out that the decay rate for the Klein-Gordon model is i^"/^ in 
contrast to the rate t~("-~'^)/2 for the free wave equation. There exist uniform 
estimates for solutions to the damped wave equation, which happen to give 
the same decay rate of i~"/^ . 

4. Are there other asymptotic representations for solutions? As example we refer 
to Nishihara [33] for an asymptotic description of solutions to the damped 
wave equation in terms of solutions of an associated heat equation and some 
exponentially decaying free waves. 

We will come back to some of these questions for more general hyperbolic models 
in the subsequent sections. An extensive study of higher order hyperbolic equations 
with constant coefficients and their dispersive estimates was done in Ruzhansky- 
Smith [49], we restrict ourselves here in this informal introductory part to the 



wave equation (1.7|. Solutions in C (M;,y'(M")) can be studied in terms of the 
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partial Fourier transform with respect to the x-variable. This gives the parameter- 
dependent ordinary differential equation 

uu + \^fu^0 (f.l2) 

with solutions explicitly given as 

u{t,0 - cosm)uo{0 + ^^^^Si(0- (1-13) 

Estimates in L^-scale can be obtained immediately from Plancherel identity. We 
will only give the following lemma stating higher order energy estimates for the 



free wave equation. The proof is straightforward from (1.13) 



Lemma 1.2. The solutions to the Cauchy problem ( 1.1 1 satisfy the a priori estimate 

\Ht.-)\\L^ <\\uQ\\L^+t\\u4L-^ (1.14) 

together with 

llD^^DXi, ■)\\l2 < Ck,a (ho||ff^+i«i + ||ui||h.+i»i-i) (1.15) 

for all fc e No and a G Ng, |a| > 1. 

The proof of the dispersive estimates of Theorem |1.1| are based on the repre- 
sentation of solutions as sum of Fourier integrals 

.i(^-«±*l«l) 



eH-^±'l4l^a,-±(t,0ii,(0dC (1.16) 

with amplitudes aj^±{t,^) supported in t|^| > 1 and an (easy to obtain) better 
estimate for the remainder term localised to i|^| < 1. A dyadic decomposition of 
frequency space is used to reduce this to integrals over spherical shells in ^ and 
they are treated by means of usually stationary phase estimates. 

1.3. Stationary phase estimates 

We will conclude this introductory section by reviewing some estimates for Fourier 
transforms of surface carried measures and relate them to estimates for Fourier 
integrals. Such estimates have a rather long history and go back to the original 
work of van der Corput [7] , [8] , Hlawka [15] , [16] , Randol [39] and in particular 
for applications in hyperbolic partial differential equations Strichartz |54j , [56] and 
Littman [55] . Later, Sugimoto [ST] , [35] , established further decay rates depending 
on the geometry of the level sets of the phase, in the analytic case, introducing 
notions of convex and non-convex indices in this context. We will rely on further 
extension of these by Ruzhansky [48], [53], allowing phases and amplitudes of 
limited regularity and depending on parameters, with uniform estimates with the 
respect to parameter, bridging the gap between the van der Corput lemma and 
the stationary phase method. 

Let S be a smooth closed hypersurface embedded in E" . We assume first that 
S encloses a convex domain of K", which is in particular the case if the Gaussian 
curvature of S is non- vanishing. Let further / e C°°(S]) be a smooth function 
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defined on the surface E. In a first step we are interested in decay properties of 
the inverse Fourier transform 

/» = / e'^-«/(Ode, (1.17) 

the integral taken with respect to the surface measure induced from the ambient 
Lebesgue measure. Before stating the resuh, we introduce the contact index 7(S) 
of the surface E to be the maximal order of contact between E and its tangent 
lines. For a 2-plane H containing the normal of the surface E at p S E we define 
7(E;p, H) to be the order of contact between the tangent TpE n H and the curve 
E n i/ in the point p. Consequently, we set 

7(E) =max max -v(T.:p,H). (1.18) 

pes H-.NpScH 

Then the following result is valid. 

Lemma 1.3. Assume E is convex. Then the estimate 

|/»|<(x)-w^||/||c'»(s) (1-19) 

is valid for all f e C'=(E), fc > ^ + 1. 

This result is closely related to the following multiplier theorem which we 
formulate in the form due to Ruzhansky |35J HHj (see also Ruzhansky-Smith [15] 
and Sugimoto iSZi)- Let x G C°°(]R") be an excision function, i.e., we assume it is 
equal to 1 for large |^| and vanishes near the origin. We then consider the operator 

T^,r:u{x)^ [ e'(-«+*(«»|^r'-x(OS(Ode (1-20) 



for a given real- valued phase function (/>{£,) S C°°(M" \ {0}) being positively ho- 
mogeneous (pipS,) = p4>{S.)^ p > 0. The L'P-L^ -boundedness of such an operator is 
related to the geometry of the Fresnel surface / level set 

E - {^ e M" : 0(0 - 1}. (1.21) 

It is sufficient to prove the limit case for p = 1, boundedness between dual Lebesgue 
spaces follows by interpolation with the (obvious) boundedness in L^(K"). First 
we recall the definition of Besov spaces. Let </) G C^(M+) define a Littlewood- 
Paley decomposition, i.e., be such that J2jez'^('^~''^') ~ ^ ^'^^ ^^^ t ^ 0. Then 
the homogeneous Besov norm of regularity s > and with Lebesgue-indices p,q E 
[l,oo] is defined by 



V 






{2n\^-'^{2-^^\)dimL.)\ (1.22) 



while the Besov space B^ qO^^) contams all functions with 



CXD 



kill.. = ll=^-V(0^(OIIL + E (2^1l^-V(2-^|el)^(OllL.)' < oo, (1.23) 
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where V(s) = 1 - E^i0(2"^'s)- It follows 5^,2(1^") = B^aC^"^) n LP{W') and 
there are continuous embeddings L^(IR") ^-> 5° 2(1^") for 1 < p < 2 and similarly 
Bp, 2 (Hi") ^-> LP (M") for 2 < p' < 00. The opposite embeddings require to pay a 
small amount of regularity (see Runst-Sickel jlj for more details) . 



Lemma 1.4. Assume that E defined by (1.21) is convex. Then T^_r is a bounded 
operator mapping S?2(K") ^ £°°(M") for all r > n - ^. 

The statement of this lemma immediately implies dispersive type estimates 
for Fourier integral representations of the form 

u{t,x) = (27r)-" /"e'("-«+*'^(«»Wo(Ode (1-24) 

Indeed, by rescaling 77 — t^, y ~ x/t, we see that 

h(i,-)llL~<t"^ho|lB. (1.25) 



ioT r — n — ^W)- The above estimate (1.25) is of interest to us because (1.24) 
solves the Cauchy problem 

uu + 4'''{^^)u = Q, m(0,-)=mo, Ut(0,-)-0 (1.26) 



generalising ( 1.1 ). Here and later on we use the notation 4)(Dx) ■ u H> ■^^^[(f>{£,)'^iO] 
for Fourier multipliers. 

In general there is no reason for the Fresnel surface of a phase to be convex. 
Of particular importance for the study of hyperbolic systems arising in crystal 
acoustics or elasticity theory or for higher order scalar equations are generalisations 
of the Lemmata |1.3| and |1.4| to non-convex surfaces. To consider this situation, we 
define the non-convex contact index 

7o(S) = max min 7(S;n, iJ) (1.27) 

pes H-.NpScH 



as maximal order of minimal contact. The replacement to Lemma 1.3 
Lemma 1.5. For general (non-convex) surfaces S the estimate 

|/»|<(x)-^5T5)||/||ci(E) (1.28) 

is valid for all f £ C^(S). 



Again we obtain a multiplier theorem for the operator T^^^ defined in (1.20). 
Lemma 1.6. The operator T^^r is bounded B?^2(IK") ^ i°°(]R") forr>n- ;^^. 



Later on in Section [4~7| we will generalise the above estimates to include phase 
functions depending on both t and a;- variables suitable for the treatment of more 
general hyperbolic equations. 
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2. Equations with constant coefficients 

Before we pursue the analysis of equations and systems with time-dependent coef- 
ficients it is instructive to understand what happens in the case of equations with 
constant coefficients. One of the very helpful observations available in this case is 
that after a Fourier transform in the spatial variable x we obtain an ordinary dif- 
ferential equation with constant coefficients which can be solved almost explicitly 
once we know its characteristics. This works well for frequencies where the char- 
acteristics are simple. If they become multiple, the representation breaks down 
and other methods are required. In the presentation of this part we follow |49j to 
which we refer for the detailed arguments and complete proofs of the material in 
this section. 

2.1. Formulation of the problem 

Let us consider the Cauchy problem for a general scalar strictly hyperbolic operator 
with constant coefficients, 

m m— 1 

J = l ;=0 \a\+r=l (2-1) 

d\u{Q,x) ^ fi{x) (^ C^{W), ; = 0,...,m-l, .tgM". 

The symbol Pj{C) of the operator Pj{T)x) is assumed to be a homogeneous poly- 
nomial of order j, and the Ca^r are (complex) constants. We denote by L{T>t, T)x) 
and Lm{^)t^^x) the operator in (|2.1[) and its principal part, respectively. 



L„(Dt,D,) = or + ^P,(D,)D 



t 1 



-1 



(2.2) 



i(D,,D,)=L™(Dt,D,)-f ^ Y, c„,,D^D;. 

1=0 |Q|-|-r=/ 



We denote by fk{£,) and Tfe(^), k = 1, . . . ,m, the characteristic roots of Lm and 
i, respectively. Hence, we write 

m m 

Lm{r,i)^X{{T~^km, L{T,O^X{{T-Tk{i)). (2.3) 

k=l k=l 

The strict hyperbolicity of L{Dti^x) means that the functions </?fc(0 are real and 
distinct, for all ^ ^ 0. It follows that (/3fc g C°°(M"\0) and that the ^pk are positively 
homogeneous of order one, (pk{\C} = ^fkiO for all A > and £, ^ 0- 

In the general analysis it is usually desirable to have conditions on the lower 



order terms for different rates of decay of solutions to (2.1 1. However, it is more 
convenient to prove the results making assumptions on the characteristic roots. 
Some results deducing properties of characteristics from properties of coefficients 
are available, and we refer to 1491 for some details in this direction. 
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It is natural to impose a stability condition excluding an exponential growth 
of solutions in time. Namely, we assume that for all £_ E M" we have that 

ImTfc(O>0 for fc = l,...,m. (2.4) 

In fact, certain microlocal decay estimates are possible even without this condition 
if the supports of the Fourier transforms of the Cauchy data are contained in the 



set where condition (2.4) holds. However, this restriction is only technical so we 
may assume (2.4) without great loss of generality since otherwise no time decay of 
solutions can be expected, which becomes clear from the following representation 
of solutions. 

2.2. Combined estimates 



By taking the Fourier transform of (2.1) with respect to x, solving the resulting 



ordinary differential equation, and taking the inverse Fourier transform, we obtain 



that the solution to the Cauchy problem (2.1) can be written in the form 

rn — 1 

7.(t,a;)=^i?,(t)/,(x), (2.5) 

where the propagators Ej{t) are defined by 

E,{t)f{x)^ e'-«(^e-'^'«)*4(i,C))7(C)rf?, (2.6) 

with suitable amphtudes A'^{t,£). More precisely, for each k and j, the functions 
A*^(t,^) are actually independent of t on the set of simple characteristics 

^fe:-{^eM":rfc(O^Ti(0 ^l^k}. (2.7) 

For these ^?(C)i ^^ have 

Lemma 2.1. For ^ E Sk we have 

(-ip E n ^.,(0 

n (r/(C)-r,(0) 

where ^ means sum over the range indicated excluding k. Furthermore, for each 
j = 0, . . . ,m — 1, and k = 1, . . . ,m, we have that the amplitude j4^'(^) is smooth 
in Sk and 

4(0 = om~') «^ 1^1 ^ °°- (2-9) 

It turns out that it is sensible to divide the considerations of how character- 
istic roots behave into two parts: their behaviour for large values of |^| and for 
bounded values of |^|. The reason is that the large time behaviour of solutions 



of (2.1) is governed by the properties of characteristics r^, of the full equation. 
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In turn, these are similar to those of tpk for large |^|. It turns out that the key 
properties of characteristics to consider are 

• the behaviour for large frequencies in which case the roots r^ must be distinct; 

• the multiplicities of roots (this only occurs for bounded frequencies); 

• whether roots lie on the real axis or are separated from it; 

• how roots meet the real axis (if they do) ; 

• the properties of the Hessian of the roots, V?rfc(^); 

• convexity-type conditions guaranteeing better decay rates. 



For frequencies away from multiplicities we can actually establish estimates 



for the corresponding oscillatory integrals that contribute to the solution in (2.6). 



Around multiplicities we need to take extra care of the structure of solutions. This 
can be done by dividing the frequencies into zones each of which will give a certain 



decay rate. Combined together they yield the total decay rate for solution to (2.1 ). 
We now state the summarised result as it appeared in [i5] . 



Theorem 2.2. Suppose u — u{t, x) is the solution of the to*'' order linear, con 
stant coefficient, strictly hyperbolic Cauchy problem (2.1 1. Denote the characteris 
tic roots of the operator by ti(^) 
k — 1, . . . , TO, and all ^ Cz i 
which take values as follows 



, Tfn (^) , and assume that Im Tk (^) > for all 
We introduce two functions, K^^'it) and K^ \t). 



I. Consider the behaviour of each characteristic root, Tfc(^), in the region \£\ > 
M , where M is large enough. The following table gives values for the function 
Kj, (t) corresponding to possible properties of Tk{^); if Tk{£,) satisfies more 
than one, then take Kj, (t) to be function that decays the slowest as i —> oo. 



Location of Tfe(^) 


Additional Property 


K^^it) 


away from real axis 




e '**, some S > 


on real axis 


detV|rfc(C)^0 

rankV|Tfc(0 ^n-1 

convexity condition 7 

no convexity condition, 70 


71/1 In 

11 - 1 / 1 In 
t 2 l-p <ji 

11 - 1 / 1 In 
t y (p q> 



KV{t). 



Then take K^^\t) = max^^i ...,„ icf ^( 
II. Consider the behaviour of the characteristic roots in the bounded region |^| < 
M; again, take K^ '{t) to be the maximum [slowest decaying) function for 
which there are roots satisfying the conditions in the following table: 
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Location of Root(s) 


Properties 


i^W(i) 


away from axis 


no multiplicities 
L roots coinciding 


e^"**, some 5 > Q 


on axis, 
no multiplicities * 


detV|Tfe(O^0 

convexity condition 7 

no convexity condition, 70 


n-l / 1 1 \ 
t 7^9^ 


on axis, 
multiplicities*,** 


L roots coincide 
on set of codimension £ 


^L-l-t 


meeting axis 
with finite order s 


L roots coincide 
on set of codimension £ 


ti-i-!(i-|) 



* These two cases of roots lying on the real axis require some additional regularity 

assumptions; see corresponding microlocal statements for details. 

** This IS the L^ — L°° rate in a shrinking region; see |49] for the details. 

Then, with K{t) = max (_R''^^(i), _R''"(t)), the following estimate holds: 

7n — 1 

\\d^dlu{t, .)|k, < C^^rK{t) J2 WMl^.^^, ' (2-10) 

where 1 <p<2,pq — p + q, and Np — Np{a, r) is a constant depending upon p, a 
and r. 

Let us now briefly explain how to understand these tables. Since the decay 
rates do depend on the behaviour of characteristic roots in different regions, in 



Theorem 2.2 we single out properties which determine the final decay rate. Since 
the same characteristic root, say Tfc, may exhibit different properties in different 
regions, we look at the corresponding rates K''^^{t),K'--^^{t) under each possible 
condition and then take the slowest one for the final answer. The value of the 
Sobolev index Np = Np{a, r) depends on the regions as well, and it can be found 
from the relevant microlocal statements that exist for each region. 

In conditions of Part I of the theorem, it can be shown by perturbation 
arguments that only three cases are possible for large ^, namely, the characteristic 
root may be uniformly separated from the real axis, it may lie on the axis, or it may 
converge to the real axis at infinity. If, for example, the root lies on the axis and, in 
addition, it satisfies the convexity condition with index 7, we get the corresponding 
decay rate K'^'^\t) = t ^^p^q'. Indices 7 and 70 in the tables are defined as 
the maximum of the corresponding indices '){^\) and 70 (S^), respectively, where 
Ea = {C : Tk{i) — A}, over all k and over all A, for which ^ lies in the corresponding 



region. The indices 7(Sa) and ^q{Y.x) are those defined in Section 1.3 



The statement in Part II is more involved since we may have multiple roots 
intersecting on rather irregular sets. The number L of coinciding roots corresponds 
to the number of roots which actually contribute to the loss of regularity. For 
example, operator {&f — Ax){df — 2A2,) would have L — 2 for both pairs of roots 
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±1^1 and ±-\/2|f|, intersecting at the origin. Meeting the axis with finite order s 
means that we have the estimate 

distil, ZkY < cllrnnm (2.11) 

for ah the intersecting roots, where Zk = {C '■ ^n^Tk{£,) = 0}. In Part II of Theo- 



rem 2.2 the condition that L roots meet the axis with finite order s on a set of 
codimension £ means that aU these estimates hold and that there is a (C^) set M 
of codimension £ such that Zi- C M for all corresponding k. 

In Part II of the theorem, condition ** is formulated in the region of the size 
decreasing with time: if we have L multiple roots which coincide on the real axis 
on a set Ai of codimension £, we have an estimate 

m — 1 

\u{t,x)\<C{l + tf-'-'y\\M\L., (2.12) 



if we cut off the Fourier transforms of the Cauchy data to the e-neighbourhood Ai"^ 
of Ai with e = l/t. Here we may relax the definition of the intersection above and 
say that if L roots coincide in a set A^ , then they coincide on a set of codimension 
£ if the measure of the e-neighbourhood A^*^ of A4 satisfies jA^*^] < Ce^ for small 
e > 0; here M" = {^€R'' : distil, M) < e}. 

We can then combine this with the remaining cases outside of this neighbour- 
hood, where it is possible to establish decay by different arguments. In particular, 
this is the case of homogeneous equations with roots intersecting at the origin. 
However, one sometimes needs to introduce special norms to handle L^-estimates 
around the multiplicities. We leave this outside the scope of this review, and refer 
to [32] for further details, as well as for microfocal improvements of some of the 
estimates under certain more refined assumptions. 

2.3. Properties of hyperbolic polynomials 

Here we collect some properties of hyperbolic polynomials. We start with the 
property of general polynomials that roots have bounds in terms of the coefficients. 

Lemma 2.3. Consider the polynomial over C with complex coefficients 

m 

z" + ciz"-i + • • • + c™_iz + c^ = Y[{z- Zk). (2.13) 

fc=i 

// there exists M > such that \cj\ < AP for each j = 1, . . . ,m, then \zk\ < 2M 
for all k = 1, . . . ,m. On the other hand, if there exists N > such that \cj\ < N 
for each j — 1, . . . ,m, then \zk\ < max{2, 2A^}. 

Now we give some properties of the characteristic roots. 

Proposition 2.4. Let L = L(Dt,I}x) be a linear m order constant coefficient 
differential operator in Df with coefficients that are pseudo-differential operators 
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in X, with symbol 

m m 

L{r, = r" + ^ PAOr-^-' + E ^dOr"-' , (2.14) 

where Pj(A^) = X^PjiO for all A > 1, \S,\ > 1, and a^ e S^-', for some e > 0. 

Then each of the characteristic roots of L, denoted ti(^), . . . ,Tm(C); *■' con- 
tinuous in K". Furthermore, for each k — I, . . . ,m, the characteristic root t^.{S.) 
is smooth away from multiplicities, and analytic if the operator L{T)t,^x) is dif- 
ferential. If operator L(Dt,Dx) is strictly hyperbolic, then there exists a constant 
M such that, if \^\ > M then the characteristic roots ti(^), . . . , Tm(^) of L are 
pairwise distinct. 



Proof. The first part of Proposition |2.4| is simple. For the second part we follow 
[i^ and use the notation and results from Chapter 12 of [TT] concerning the 
discriminant Ap of the polynomial p{x) = Pmx"^ + • ■ ■ + Pix -\- po, 

Ap EE A{po, ...,Pm):= (-1)'"'^^^^"-' llix^ - x,f , (2.15) 

where the Xj (j = 1, . . . , m) are the roots oi p{x). We note that Ap is a continuous 
function of the coefficients Pq, . . . ,Pm of p(a;) and it is a homogeneous function of 
degree 2m — 2 in them. In addition, it satisfies the quasi-homogeneity property: 

A(po, Api, A2p2, . . . , A"V,„) = A"(™-i)A(po, . . . .Pm)- (2.16) 

Clearly, Ap = if and only if p{x) has a double root. We now write L{t, ^) in the 
form 

i(r, i) = i™(T, i) + ai(e)r"-i + aalOr""" + • ■ ■ + a,„-i(Or + a„.(0, (2.17) 
where 

m 

i„(r,e) = r™ + E^.(OT"-^' (2.18) 

is the principal part of L(r, ^). Note that the Pj (^) are homogeneous of degree j and 
the aj(^) are symbols of degree < j. By the homogeneity and quasi-homogeneity 
properties of A^, wc have, for A 7^ 0, 

Ai(AO = A(P„,(AO + a,n(AO, . ■ ■ .Pii^i) + ai(Ae), 1) 
= A(A™[P™(e) + ^^#], . . . , A[Pi(0 + ^], 1) 

^ y„(2m-2)^(p^(^) _^ a,^ ^ ^ _ . , A~(™-1) [Pi (C) + ^^^], A"™) 

(using that A is homogenous of degree 2m — 2) 

= A'"("-i)A(P,„(0 + ^-^, . . . ,Pi(e) + ^, 1) 

(by quasi-homogeneity) . 
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Now, since L is strictly hyperbolic, the characteristic roots Lpi{S,), ■ ■ ■ , fmiO of Lm 
are pairwise distinct for ^ 7^ 0, so 

AlJO - A(P„(0, . . . , Pi(0, 1) ^ for ^ ^ 0. (2.19) 

Since the discriminant is continuous in each argument, there exists 6 > such that 
if 1^47^1 < (5 for ah j = 1, . . . , m then 

|A(P™(0 + ^i^ , . . . , Pi(0 + ^ , 1)1 ^ 0, (2.20) 

and hence the roots of the associated polynomial are pairwise distinct. Since qj are 
symbols of order < j, it follows that for |^| = 1 and large A, we have Ai(A^) 7^ 0, 
which means that the characteristic roots of L are pairwise distinct for |A^| > M 
for some constant M > 0. D 

Therefore, we see that for large frequencies, the characteristic roots of the 
full symbol are smooth. In fact, for large frequencies again, we can regard the 
lower order terms of the operators as a perturbation of its principal part. Using 
this, we can show that the roots actually have a good symbolic behaviour for large 
frequencies with a certain dependence of the order of the symbol on the order of 
the lower order terms. 

We will use the standard notation for the symbol class S"^ of all amplitudes 
a = a(x,0 e C°°(M" x M") satisfying 

|af9|'a(x,C)|<C„^(0^-|"l, (2.21) 

for all multi-indices a,/3 and all x,^ G M". Here, as usual, (^) = Vl+TCF- 

Proposition 2.5. Let L = L{Dt,^x) be a hyperbolic operator of the form 

m m 

L(D„D,)=DJ"+^P,(D,)Dr'+^ J2 c„,,(D,)Dr^ (2.22) 

j=l j = l \a\+rn-j=K 

where Pj{XO = >^^ PjiO for A > 1, |^| > 1, and c„j e 5l"l. Here 0<K<m-l 
is the maximum order of the lower order terms of L. Let ti(^), . . . , TVn(^) denote 
its characteristic roots. Then there exists a constant C > such that 

|Tfe(C)|<C(0 forall(,eR",k = l,...,m. (2.23) 

Suppose in addition that L is strictly hyperbolic, and denote the roots of the prin- 
cipal part Lm{T, ^) by ipi{S^), . . . , (^m(^). Then the following holds: 

1. For each t^^S,), k — 1, . . . , to, there exists a corresponding root of the principal 
symbol (fkiO {possibly after reordering) such that 

kfe(0-^fe(OI<C(e)^+i-" for all ^eR". (2.24) 

2. There exists M > such that, for each characteristic root of L and for each 
multi-index a, we can find constants C — Ck.a > such that 

151^x^(01 < C|er-I"l forall\i\>M, (2.25) 
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3. There exists M > such that, for each Tfc(^) a corresponding root of the 
principal symbol ipk{S,) can be found {possibly after reordering) which satisfies, 
for each multi-index a and k — 1, . . . ,m, 

\d?rk{0 - a.>fc(C)l < C|e|^+^-"-l"l for all |^| > M (2.26) 



The estimate (2.23) follows immediately from Lemma |2.3| The symbolic be- 
haviour can be shown by perturbative arguments. In particular, since we always 
have X < m — 1, we get the following special case of Proposition |2.5| which we 
formulate for its own sake: 



Corollary 2.6. Let L{Dt, Dj.) be a strictly hyperbolic operator as in Proposition 2.5 
Then we have 

\rk{0-MO\<C for all ^eR", (2.27) 

and there exists constants M > and C > such that we have 

|Vt,(0|<C forall\S,\>M, (2.28) 

and 

\d^rk{0 - d1M0\ < C|er'"' for all \^\ > M , (2.29) 

for each multi-index a and fc = 1, . . . , m. 

2.4. Estimates for oscillatory integrals 

Theorem |2.2| is composed of a number of statements each giving a corresponding 
time decay rate in a corresponding frequency region. For the precise formula- 
tions and proofs of these statements we refer to j49]. However, anticipating the 
estimates for the oscillatory integrals appearing in the analysis of systems with 
time-dependent coefficients in subsequent sections, we give now two estimates cor- 
responding to large frequencies in which case the properties of hyperbolic polyno- 
mials in Section [2.31 can be used. 

We start with the case where all the level sets 

S^(t) = U G M" : t{0 = A} (2.30) 

of the phase function r : M" — )■ M, if non-empty, are convex in the sense that the 
set enclosed by Ea is convex. For simplicity of the formulation, we agree that the 
empty set is also convex. 



We also note that if t is a characteristic of the operator in (2.1), we have 



to estimate the oscillatory integrals appearing in (2.6). These are dealt with in 



Theorem |2.7| below. Proposition |2.5| assures us that the function t satisfies the 



assumptions of Theorem 2.7 for large frequencies ^ 



Theorem 2.7. Let r : M" — > M &e such that Ea(t) are convex for all A > 0. 
Let X S C°°(M"), and assume that on its support suppx we have the following 
properties: 

1. for all multi-indices a there exist constants Ca > such that 

|a|'r(0|<C„(0^-l"l; (2.31) 
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2. there exist constants M, C > such that for all |^| > M we have |t(^)| > C|^|; 

3. there exists a constant C > such that \d^T[\u])\ > C for all ui G §"^^, 
A > 0; in particular, |Vt(^)| > C for all ^ € R" \ {0}; 

4. there exists a constant i?i > such that, for all A > 0, 

X-'j:xiT) C Br,{0) . (2.32) 

Set 7 :~ sup_)j^Q 7(S^(t)) and assume this is finite. Let Oj = aj{£^) G S^^ he a 
symbol of order — j on M". Then for all t > we have the estimate 



J(x<+r(i)t) 



a.iOxiOfiOd^ 



<c(i + f)-'^(^-i)||/IUp 



(2.33) 



Li 



where pq = p + q, 1 < p < 2, and the Sobolev order satisfies Npj^t ^ 't-(- ~ ~) ^ J 
forO<t<l,andNp^^,t>{n~^)il-l)-jfort>l. 



Estimate (2.33) follows by interpolation from the L^ — L^ estimate (easy part) 



combined with an L^ — L°° estimate (harder part). The estimate for large times 



follows from the L°°-estimate for the kernel of the integral operator in (2.33). In 
fact, for the kernel we have the following 



Proposition 2.8. Under the conditions of Theorem 2. 7 with x = 1 and with the 
assumption that a G C^(M") the estimate 



follows for all x £ 



A^<+T{m 



and all < > 0. 



a{Od^ 



<c{i + ty— 



(2.34) 



We now give an analogue of Theorem |2.7| without the convexity assumption 
on the level sets of t. In this case we can only assure a weaker, one-dimensional 
decay rate. 

Theorem 2.9. Let r : M" — > M 6e a smooth function. Let x G C°°(M"), and on its 
support supp X; assume the following: 

1. for all multi-indices a there exist constants Cq > such that 

19^(01 <C„ (0^-1"!; (2.35) 

2. there exist constants M, C > such that for all |^| > M we have |t(^)| > C|^|; 

3. there exists a constant C > such that |9^T(Aaj)| > C for all ui G §"~^ and 
A>0; 

4. there exists a constant Ri > such that for all X > we have 



x-'j:x{t) c Br,{o) . 



(2.36) 



Set 7o := sup_^^q7o(Sa(''")) and assume it is finite. Let aj = aj(^) Cz S ^ be a 
symbol of order —j on R". Then for all t > we have the estimate 



3'(-«+-(«)*)a,(0x(0/(0'^e 



<c(i + <)-*('-')||/ILp 



(2.37) 



Li 
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where pq ^ p + q, 1 < p < 2, and the Sobolev order satisfies Npj^t ^ ?^( ) ^ J 

forO<t<l, and Np,,,t > {n - ^)(i - i) - j for t > 1. 

As in the convex case, the main thing is to show the fohowing kernel estimate. 

Proposition 2.10. Under the conditions of Theorem \2.S\ with x = 1 o,nd the as- 
sumption that a G C^(M") the estimate 






<C(l + i)"^ (2.38) 



holds true for all x £ M" and t>Q. 

As a corollary and an example of these theorems, we get the following possi- 
bilities of decay for parts of solutions with characteristic roots on the axis. We can 
use a cut-off function x to microlocalise around points with different qualitative 
behaviour (hence we also do not have to worry about Sobolev orders). Combining 
the above theorems with what we can obtain by the stationary phase method, we 
get 

Corollary 2.11. Let J7 C M" he an open set and let t : H. ^ M. be a smooth 
real-valued function. Let x £ C^{Q). Let us make the following choices of K{t), 
depending on which of the following conditions are satisfied on suppx- 

(1) //detVV(0 ^0 for all ^ e Q, we set K{t) = (l + t)"?^^"^^ 

(2) //rankVV(^) = n - I for all £, e n, we set K(t) = {l + ty'^^p^^\ 

(3) // T satisfies the convexity condition with index 7, we set K{t) = (1 -|- 

t) 'i ^ P 1 \ 

(4) If T does not satisfy the convexity condition but has non-convex index 70, we 
set K{t) = {l + ty^^p~^\ 

Assume in each case that other assumptions of the corresponding Theorems |2.7|- 



2.9 in cases (3)-(4) are satisfied, as well as assumptions in |49j for cases (l)-(2). 



Let 1 < p < 2, pq ^ p + q. Then for all t > we have 

ei(-?+-(e)*)„(^)^(^)7(^)d^ <CK{t)\\f\\LP. (2.39) 

Li 

We note that no derivatives appear in the L^-norm of / because the sup- 
port of X is bounded. In general, there are different ways to ensure the convexity 
condition for r, some of the criteria given in [49) . 

3. Some interesting model cases 

We now look at equations with time-dependent coefficients. In this section we will 
review two scale invariant model cases, which can be treated by means of special 
functions. They both highlight a structural change in the behaviour of solutions 
when lower order terms become effective. This change is a true variable coefficient 
phenomenon, it can not arise for equations with constant coefficients as treated in 
Section O 
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3.1. Scale invariant weak dissipation 

We follow the treatment in :63j and consider the wave model 

wtt - Am+ — Ui == 0, u(l,-)=uo, Uf(l,-)=Mi (3.1) 

for data uq, ui £ ^'(M"). By partial Fourier transform the above problem reduces 
to the ordinary differential equation 

utt + \^\^u+'^ut^0 (3.2) 

parameterised by |^p. We will proceed in two steps, first we construct a fundamen- 
tal system of solutions to this equation and after this use it to provide a Fourier 



multiplier representation to solutions of (3.11 



3.1.1. Reduction to special functions. Equation (3.2) can be reduced to Bessel's 



differential equation. Let p = i — /x. If we look for particular solutions of the form 



^(^:C) — (i|^l)''^(^l^l)i a short calculation reduces (3.2) to 

T^v" + Tv' + (r^ - p^)v = 0, (3.3) 

which is Bessel's differential equation of order p. For large values of r we will use 
the fundamental system of solutions given by the Hankel functions H^ . In order to 
get a precise description of its solutions for small r, we have to distinguish between 
integral p and non-integral p. In the first case, solutions can be represented by the 
fundamental system J7p, yp, while in the latter we use J^±p. 



p^Z, 


(3.5) 


peZ, 


(3.6) 



Lemma 3.1. Any solution to ( |3.2[ ) can be represented in the form 

uit,o^c+iomrKit\^\)+^-(omrK(i\^\)' 
= Aiomirj.pm) + Bmmrjpmi 
= A{omrj-pm\) + mmry^pm), 

with coefficients depending linearly on the Cauchy data 

c±{0 = Ci^oiOMO + c±AOMO (3.7) 

and similarly for A{^) and B{^). 

ExpHcit expressions for C±_j{^), Aj{^) and Bj{^) can be obtained from these 
formulas and known expressions for derivatives and Wronskians of Bessel functions. 
Wc will not go into these details here and refer to the reader to the original paper 
[55] . Of greater importance for us are consequences about structural properties of 
solutions. They are based on elementary properties of the Bessel functions. We 
collect some of them in the following proposition. They are taken from the treatise 
of Watson, (BUI Sections 3.52, 10.6 and 7.2]. 

Proposition 3.2. 1. The functions "H^ (r) possess an asymptotic expansion 

oo 
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Figure 1 . Decomposition of the phase space into zones 



as r —> oo which can be differentiated term by term; 

2. the function t^''J^p{t) is entire and non-vanishing in t — 0; 

3. we have 

2 

yn(T) = -Jn{T)l0gT + An{T) 

with T^-AniT) entire and non-vanishing in t — 0. 



(3.9) 



If we want to derive properties of solutions, we have to distinguish between 
large i|^| and small i|^|. This distinction will later on lead to the introduction of 
zones and will also play a crucial role in the definition of symbol classes and in the 
more general consideration of hyperbolic systems. 

3.1.2. High frequency asymptotics. The asymptotic expansion for Hankel func- 
tions immediately yield that 

?^(t,e) ^C±(Oe±'*l«l (a±t-ner^ + l-o.t.) (3.10) 

as t|^| — > oo. A more precise look at the terms C'±{^) gives a representation 

C±A0-\Cr'+^-O-t., .?=0,1, (3.11) 

such that if |^| > 1 and t\^\ > 1 both can be combined to the following (very rough) 
description. The appearing Fourier multipliers have a uniform decay rate i~^ which 
corresponds to the high frequency energy estimate of the following lemma. 
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Lemma 3.3. Assume ^ suppuo andO ^ suppwi. Then the solutions to the weakly 



damped wave equation (3.1) satisfy 

Mt,x)\\L2 <t"'' iWuoh^ + \\u4h-^) (3.12) 

and the energy estimate 

E{u;t)<t-^^'E{u;0). (3.13) 

The constants in these estimates do in general depend on the distance of 
to the Fourier support of the data. 

We can be shghtly more precise than the above lemma. If we muhiply u{t, S.) 
by f^ the main term reduces to a representation of a free wave while lower order 
terms decay at least as t^^. This gives a description of the large-time asymptotic 
behaviour of weakly damped waves in terms of free waves, see also [64j. 

Lemma 3.4. Assume ^ supp Uq and ^ supp Ui . Then there exists a free wave 
w, i.e., a solution to the free wave equation 

Wft - Aw = 0, w{l,-)^wo, wt{l,-)^wi (3.14) 



to appropriate data, such that solutions to the weakly damped wave equation (3.1) 
satisfy 

\\t^'u{t, •) - w{t, ■)\\l2 < t-' (WuoWl- + \\u,\\h-^) (3.15) 

The operator assigning the data wq, wi to uq and ui is bounded on L^ x H^^ . 

Proof. The proof is based on the elementary expression of 'w{t, ^) in terms of the 
inital data, 

2€}it,0 = e"l«l(zBo(C) " l^r''^i(C)) +e-"l«l(^o(0 + l^r'^i(O)- (3.16) 
Therefore, we relate the initial data wq and wi to uq and ui by the system of 
linear equations 

MO T m-'MO = 2|^r^(c±,o(OSo(C) + C±^i{OMO)- (3.17) 

When forming now the difference (|3.15[) this cancels the main terms of (3.101 and 



we are left with terms decaying like t ^ or faster. Note that (3.17) determines the 



data Wq and wi and yields the desired boundedness property. D 

3.1.3. Low frequency asymptotics. The situation for low frequencies is completely 
different. We restrict ourselves to the case of non-integral p (and except for /i = 1 
this will not alter any estimates we provide here). Solutions are represented by 



(3.5). Crucial point is that we get no decay in time for the multiplier (i|^|)''j7_p(i|^|), 
while the behaviour of (i|CI)''»//9(^ICI) depends on whether p > or p < 0, i.e., 
whether fi < ^ or /i > |. In the first case, estimates of solutions are in general 



increasing in time (similar to (1.14)). In the second case or when estimating higher 
derivatives of solutions (i|^|)''J'_p(i|^|) will become the dominant part. Then any 
form of decay in time has to come from the behaviour of A{^) near ^ = 0. If A(^) 
vanishes to order k, we get an estimate by t~'^ uniform in i|^| < 1 and if k is not 
too large also uniform in |^| < 1. This can be used to deduce the following kind of 



higher order energy estimates for solutions to (3.1 ) 
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Lemma 3.5. Let /i > 5 ^^"^ A; G N satisfy k < fj,. Then the following higher order 
energy estimate 

WMt.■)\\L^<t-''[\\u4H^+\\ul\\H^-^) \a\ = K (3.18) 



is satisfied by any solution of (3.1). 



Sketch of proof. We will omit some of the details here, for the full argument we 
refer to [Ml Sect. 3.1]. As mentioned above, the decay rate for high frequencies 
is t^^, while the decay for small frequencies has to be related to a certain zero- 
behaviour of the coefficient A{£^) in ^ = 0. This is done by taking ij'^-norms on the 
right-hand side and estimating just homogeneous TJ'^-norms on the left-hand side. 
This gives an additional factor of |^|'' for small |^| and in turn a uniform estimate 



hy t '' within the zone i|^| < 1. Note, that the second term in (3.5) decays faster, 
as there is an additional i^~^^ factor. 

It remains to discuss the intermediate part where t\^\ is large and |^| remains 
bounded. Here we will estimate the multiplier by |CI'^(i|CI)~^ ^ ^"'^i which is true 
whenever k < ^. D 

We remark that the proof from [53" also yields a similar statement involving 
both t and x-derivatives of solutions. We decided to omit this, because then more 
detailed arguments involving recursing formulae for Bessel functions would have 
been needed here. 



3.1.4. Notions of sharpness. There are different ways to measure the sharpness 
of a priori estimates. We will mention two of them and explain their importance 
related to the further considerations. If we are given an energy estimate of the 
form 

E{u;t)<f{t) (3.19) 

it should always be equipped with a class of data. A first question should be: Can 
we find data from this class, such that the energy really does behave in this way? 
This is clearly the case of the high frequency energy estimate of Lemma |3.3[ In 
fact, for all data with frequency support away from zero this happens. If /i < 1 
then this happens even for all data from H^ x L^, see |M| . 

On the other hand, if we cannot find such data we may ask whether we can 
improve the estimate to E{u;t) < g{t) for some g{t) = o{f{t)). If this can be 
done, the estimate clearly was not sharp. However, if we do not find data with the 
prescribed rate but also can not improve the rate, a different kind of sharpness 
appears. This is seen in the low frequency asymptotics. For any /i > 1 the estimate 

||Vu(t, ■)\\l2 + \\ut{t, ■)\\l2 < t-' {\\uo\\h^ + \\ui\\l2) (3.20) 

is of this form. 



22 M. Ruzhansky and J. Wirth 

3.2. Scale invariant mass terms 

This situation was studied recently by Del Santo-Kinoshita-Reissig [9] and in the 
PhD thesis of C. Bohme, |5, following along similar lines to the treatment in the 
previous section. We will only sketch the major differences and the conclusions to 
be drawn from them. 

The model under consideration is the Cauchy problem 

utt- Au+ —u = 0, u{1,-)=uq, Uf{l,-)^ui (3.21) 

for a Klein-Gordon equation with time- dependent mass. Again we assume data 
to belong to ^'(M") and we reduce the problem by a partial Fourier transform to 
the ordinary differential equation 

^2 

utt + m^u + —u = 0. (3.22) 

This differential equation can be related to Kummer's confluent hypergeometric 
equation. Let 2/3=1 + y/1 — k^. If we look for particular solutions of the special 
form 

uit,o = e''\^\mrvm) (3.23) 

we obtain with the substitution r — 2it\^\ 

Tv" + (2p - t)v' + T-y = 0. (3.24) 

Solutions to this equation are given by confluent hypergeometric functions. A 
system of linearly independent solutions is given in terms of Kummer's functions 

by 

eo(p,2p;2ii|^|), eo(l-p,2-2p;-2ii|e|) (3.25) 

for 

eo(a,/3;r) = ^;'^'^^^;' ^^^' (3.26) 

In combination with known statements about these functions we can again describe 
the Fourier multipliers appearing in the representation of solutions. They behave 
different in different zones of the phase space, the decomposition is the one depicted 



in Figure 3.1.1 Furthermore, we have to distinguish between small values of k and 



large values. 

The case k < 1. In this case the parameter p £ [1/2, 1) is real and a solutions to 



(3.22) show the following behaviour. If t|^| < 1, the main terms of the fundamental 
solution behave like f and t^~P (and with improvement of one order for any t- 
derivative or multiplication by |^|), except for k = 1 and p — 1/2 where an 
additional log-term appears. The large-frequency behaviour on t|^| > 1 is described 
again by the oscillatory terms exp(±ii^) similar to the high-frequency expansion 



(3.10) for the scale invariant weak dissipation case. 
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The case k > 1. Now the parameter p is complex with real part equal to 1/2. 
When i|^| < 1, solutions to (3.22[) grow like t^/^ but also exhibit an oscillatory 



behaviour like exp(±i-\/K^ — Hogi). For t\S} > 1 the solutions are bounded and 
oscillating. The high-frequency expansion will have a changed phase. 



4. Time-dependent hyperbolic systems 

In this section we will provide a diagonalisation based approach to obtain the 
high-frequency asymptotic properties of the representation of solutions for more 
general uniformly strictly hyperbolic systems. The exposition is based on ideas 
from the authors' paper "ST]. 

One immediate complication compared to the constant coefhcients case in 
Section [2] is that there is no simple formula similar to (2.8) for amplitudes of the 



solutions. Thus, the first task is to construct suitable substitutes for (2.8), which 
can be done by different methods depending on the properties of the coefficients 
of the equation. In the presence of the lower order terms of the equation, an 
hierarchy has to be introduced in order to fall into the construction scheme for 
the amplitudes. Thus, in Section[2]any lower order terms could have been allowed, 
but here, we need to impose decay conditions, for them to fall into the required 
symbolic hierarchy. 

4.1. Motivating examples 

The motivation to consider systems in this framework is two-fold. On the one 
hand, the treatment of wave models and more general higher order hyperbolic 
equations naturally leads to a reformulation as pseudo-differential hyperbolic sys- 
tems. This was key ingredient for Reissig-Smith j46j to treat wave equations with 
bounded time- dependent propagation speed. In [65] , [66] , Wirth considered time- 
dependent dissipation terms and discussed their influence on energy and disper- 
sive type estimates for their solutions. The treatment of the non-effective case fits 
to the considerations presented here. Equations with homogeneous symbols and 
time-dependent coefficients appear naturally also in the analysis of the Kirchhoff 
equations. Higher order equations of Kirchhoff type were discussed by Matsuyama 
and Ruzhansky in '29] . 

On the other hand, systems are of interest on their own. Recently, D'Abbico- 
Lucente-Taglialatela [1] studied general time-dependent strictly hyperbolic differ- 
ential systems of the form 

n 

DtC/ = ^A,(t)D,^.(7 + i?(i)t/, f/(0,-) = C/o, (4.1) 

i=i 

where as usual D = —id and the matrices Aj{t),B{t) € C'*^'' satisfy natural 
symbolic and structural assumptions. Besides guaranteeing a form of uniform strict 
hyperbolicity up to t = cx) it is of importance to control the amount of oscillations 
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in coefficients. In this setting this is usually done by assuming 

Aj(t)er.{0}, B{t)e%{l}, (4.2) 

where 

r.m = 1/ e c^{^+) ■■ |Df/(t)| < Ck (^^^^TT^)'^ f ^^-^^ 

and I' g [0,1]. We will concentrate on the case v — here; this simplifies the 
consideration without omitting most of the main ideasr] Key assumption will al- 
ways be that the characteristic roots, i.e. the solutions Xj{t, ^), j = 1, . . . , d of the 
polynomial equation 

detlxit,O~J2MmA=0 (4.4) 

are real and uniformly separated from each other for t > and |^| = 1. 

The approach also applies to higher order scalar hyperbolic equations with 
purely time-dependent coefficients of the form 

m— 1 

^Tu+Y. E aj,a(i)D^D> = 0, D^w(0,-)=wj, j=0,l,...,m-l, 

(4.5) 
where 

aj,a{t)(^T,{m-j-\a\}, (4.6) 

provided the roots of the homogeneous polynomial 

rn— 1 

^'" + E E a,,o(i)A^r = o (4.7) 

j=0 |Q|=m-j 

are real and again uniformly separated for i > and |^| = 1. 

Note, that our approach uses a high regularity of coefficients (resulting in 
sharp conditions on the decay behaviour) . There is an alternative approach due to 
Matsuyama-Ruzhansky [30j based on an asymptotic integration argument applica- 



ble to homogeneous higher order equations of the form (4.51 with a' ^ (t) E L^ (IR+) 
for j -|- |a| = TO and aj,a{t) = otherwise. This yields uniform bounds on the en- 
ergy of solutions as well as it allows for the derivation of dispersive estimates. Main 
difference is that one has to use a higher regularity of the data to compensate the 
low regularity of coefficients in order to obtain sharp dispersive decay rates, see 
Section Ol 



'^We omit the index u ii u = 0. 
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4.2. Symbol classes 

The following calculations are based on a decomposition of the phase space into 
different zones. They correspond to the distinction between small and large fre- 
quencies used in the previous section, see Figure |3.1.1| We denote 

2hyp = {(i,0:(l+i)lel>l} (4.8) 

as the hyperbolic zone and set Zpd — (K+ x M") \ -Zhyp. This decomposition was 
inspired by the treatment of Yagdjian, j69j . We denote the common boundary 
between both regions by t^ and use the notation Xhyp{t,0 = x((l + ^)IC|) for a 
smooth excision function supported inside -Zhyp and being equal to one for (1 + 

m\ > 2. 

We denote by S {mi, 1112} the set of all time-dependent Fourier multipliers 

a(f, ^) e C°°(]R+ X M") satisfying the symbol estimates 



|DfDfa(t,e)| <C,,„ (^maxjiei,^ 



"ll — |ct| / 1 \ '"2+fe 



1+i 



(4.9) 



for all k and multi- indices a. The set S {mi, 1712} possesses a natural Frechet space 
structure. Furthermore, the classes behave well under forming products, taking 
derivatives and certain integrations. The excision function Xhyp belongs to 5{0, 0}. 

Proposition 4.1. 1. S{mi — k,m2 + k'} C S{mi,m2} whenever k' > k; 

2. 5{mi,r7i2} • S{m[, m.2} C S{m,i + m.'i,m,2 + m,^}; 

3. D^D^5{TOi,m2} C5{mi- |a|,m2-hfc}. 

4. If a{t,S,) e 5{— 1,2} satisfies suppa(i, ^) C 2^hyp7 then 

/CO 
a(T,Odr, Cy^O (4.10) 

satisfies Xhyp{t,i)b{t,S,) £ 5{-l, 1}. 

Proof. We only explain how to prove the last statement, the others arc obvious. 
From the symbol estimates we get for t>t^ 

and the uniform bound by C for t < t^. As derivatives with respect to t clearly 
satisfy the right estimates, we are left with ^-derivatives. For them we obtain in a 
similar fashion 



Pf'(«)l<c-f ^j^^2__a. = - 



c 



t)l^ll+kl 



whenever t > t^. Combining this with support properties of Xhyp(^jO s-i^d its 
symbol estimates the desired statement follows. D 
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Of particular importance for us will be the embedding S{mi — k, 1112 + A;} '->■ 
5{r7ii, 7712}- It is customary to denote the residual class of this hierarchy by 

n{l>}^ fl 5{mi,m2}. (4.12) 

mi+7n2— ^ 

We will also need a certain homogeneous version of these classes allowing for 
singluraties in ^ = 0. We write 5*{mi,7n2} to denote the class of all functions 

a(i,^) e C°°[R+ X (M" \ {0})) satisfying 

|DfD^"a(i,OI <Cfc.„|Cr-l"l -— (4.13) 



uniform in (1 + t)|^| > Cq for some constant Cq > 0. 

4.3. Uniformly strictly hyperbolic systems 

We collect our main assumptions now. We use the notation 5{mi,rn2} also for 
matrix-valued multipliers and denote operators corresponding to such multipliers 
by a(i,Da;), i.e., we define 

a(f,D,)/-^-i[a(t,e)7(0]- (4.14) 

Then we consider the Cauchy problem 

Dft/^ A(i,D^)C/, [/(0,-) = t/oe^'(M";C'^) (4.15) 

for a. d X d matrix A{t,£^) e 5{1,0} satisfying suitable conditions. 

(Al): There exists a positively ^-homogeneous matrix function yli(i,^) € 5*{1,0} 
satisfying Ai{t,p^) — pAi{t,S_), p > 0, such that 

Ait,0-Ai{t,OeS,{0,l}- (4.16) 

(A2): The eigenvalues Xj{t,^) of Ai{t,S^) are real and uniformly distinct in the 
sense that 

inf \'^i^^0-'^(^^0\ > (4.17) 



for all i y^ j. 



We call (4.15) uniformly strictly hyperbolic if both of these assumptions are satis- 
fied. The symbolic calculus of Proposition |4.1| allows to draw several conclusions. 
Statements and proofs are taken from [51j . 

Lemma 4.2. 1. The eigenvalues Xj{t,^) of Ai{t,S.) satisfy the symbol estimates 

X,it,Oe 541,0} (4.18) 

together with 

(A,(t,0~A,(t,0)"'e 54-1,0}. (4.19) 

2. The spectral projection Pj{t,£^) associated to the eigenvalue Xj(t,£^) satisfies 

P,(t,Oe 540,0}. 
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3. There exists an invertible matrix M{t,£^) G 5* {0,0} with inverse satisfying 
M-'^{t,i) e 5* {0,0} such that 

A^{t, OM{t, = Af (t, OV{t, (4.20) 

holds true for I?(i,^) = diag(Ai(i,^), . . . , Ad(t,0) the diagonal matrix with 
entries Xj{t, £,)■ 

Proof. 1. The properties of the characteristic roots follow from the spectral esti- 
mate |Aj(i,^)| < ||^i(i,^)|| together with the symbol properties of the coefficients 
of the characteristic polynomial and the uniform strict hyperbolicity. Indeed, dif- 
ferentiating the characteristic polynomial 

d 

= ^/rf_fc(t,e)(A,(i,O)', h{t,O^S4k,0} (4.21) 

fc=0 

yields the linear equation 

d d 

DiA,(t,C)^fc/d-.(t,e)(A,(t,0)'"' = -E(Dt^'i-fe(^'?))(^j(^'0)' (4.22) 



, fe-i 

fc=l fc=0 



for the i-derivatives. The assumption of uniform strict hyperbolicity (A2) is equiv- 
alent to a uniform bound for the inverse of the sum on the left-hand side, indeed 
it follows from Vieta's formula that 
d 

Y,kId-kit,0{^iit,Oy~' =l[{XAt^O - ^dt,0) (4.23) 

fe=i i^j 

holds true for any fixed j. This yields the desired estimate for DtXj{t,£,). Iterat- 
ing this procedure gives expressions for higher time-derivatives in terms of lower 
ones and implies corresponding estimates. On the other hand, smoothness alone 
together with the homogeneity implies the symbolic estimates with respect to ^. 
Finally, the symbolic estimates for the inverse of the difference follow from 

differentating the identity (A,(i,0 - Xjit^O)'^ {Mt,0 - Aj(i,0) = 1- 

2. The eigenprojections can be expressed in terms of the characteristic roots 
as 

such that the symbol estimates for Pj{t,£^) follow directly from the symbolic cal- 
culus of Proposition |4.1[ 

3. A symmetriser of the matrix Ai{t,^), i.e., a matrix H{t,^) such that 
H{t,£^)Ai{t,£^) is self-adjoint, is given by 

d 
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and therefore satisfies H{t, ^) e 5*{0, 0}. We can express the inverse of the diago- 
naUser in terms of this symmetriser. Let Vj{t,^) be a (smoothly chosen) unit vec- 
tor from the one-dimensional jth eigenspace range Pj(i,^) and M~^(i,^)C for any 
fixed C e C* be the vector with the inner products {vj,Pj{t,^)C,) — (i/j,H{t,^)(^) 
as entries. 

Since Vj(t,^) is unique up to sign locally in t and ^, it is expressible as 
-Pj(^j'^)C/ll^j(^jC)CII for some fixed vector C, £ C^ chosen away from the com- 
plement of the eigenspace. Then differentiating this expression directly implies 
i^jit,^) e 540,0}. This implies M-^{t,() € 5*10,0} by the above definition of 
M~^{t,^). Furthermore, A/(t,^) has the vectors Vj{t,£) as columns and again 
A'/(t,0 e 540,0} follows. D 

4.4. Diagonalisation 



As the next step we want to diagonalise the system (4.15) within the 5{-, •} hier- 



archy modulo remainders from %{!} and terms supported within Z-p^. 

AAA. Initial step. We first diagonalise the main part Ai{t,£) using the matrix 
family M{t,i) e 5*{0,0}. We denote 

V^"\t,O^M-\t,OU{t,0, (4.26) 

within 2^hyp, so that a short calculation yields 

D,y W = (I?(t, -I- Roit, 0) V^°^ (4.27) 

with 

+ (D,M-i(t,a)Af(t,0 
€540,1}. 



We will use (4.27) as starting point for a further (sequence) of transformations 



applied to the system. 

4.4.2. The diagonalisation hierarchy. We will construct matrices 

N^''\t,0 e S4-k,k} and F^'''>{t,0 e S^-k^k + l} 
such that with 

K K-1 

7VK(i,e) = i+E^^'^(*'^)' FK-i{t,o^Y.^^'''i*'0 (4.28) 

fe=l fc=0 

the operator identity 

-NKit,0{^t~T^{t,0-FK-i{t,0)^S4-K,K + l} (4.29) 
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is valid. The construction is done recursively. The identity (4.291 will yield com- 
mutator equations for the matrices N'-''\t,S,) and the matrices F('^'(i,^) are de- 
termined by their solvability condition. The strict hyperbolicity of the system is 
crucial for this construction. 



First step, K = 1. We collect all terms from (4.29) which do not a priori 



belong to 5*{— 1, 2} and require their sum to vanish. This yields the condition 

[V{t,0,N('\t,0]=-R,{t,0 + F^°Ht,0- (4.30) 

Because 'D{t, ^) is diagonal, the diagonal entries of the commutator vanish and we 
must have 

FW(i,0 = diagi?o(i,0 (4.31) 

for solvability. Furthermore, the strict hyperbolicity assumption (A2) implies that 
the solution is unique up to diagonal matrices and the entries of the matrix 
N'^^^t,^) must be given by 

The diagonal entries will be set as {N'^^'>{t, ^))^ ^ = 0- ^ is evident that i^(°)(i, & 
540, 1}, while 7V(i)(t,C) e 54-1, 1} follows from Lemma 



4.2 



Recursion fc 1— >■ fc -f 1. If we assume that all terms up to order k are already 
constructed and satisfy the symbolic inequalities then the conditions for the next 
terms follow from the requirement 

Bk+iit,0 - Bk{t,0 e S4-k - l,k + 2}, (4.33) 



with Bk as in (4.29). Indeed, collecting again just the terms which do not a priori 
belong to this symbol class and setting their sum to zero yields the commutator 
equation 

[V{t,^),N<^''+'\t,0] = B,it,0 + F^''\t,0- (4.34) 

The solution is again given by 

<"'"■'<«»« -Arofc)' ■^^- <"^' 

provided we set 

i^('=)(t,0 = -cliagBfc(t,0. (4.36) 

The diagonal entries will again be fixed as (A^*^''+^)(i,^)) . . = 0. Clearly, the di- 
agonal terms satisfy F^'''>{t,^) £ 5*{— fc, fc + 1} as consequence of the assump- 
tion Bkit,^) e 5*{-fc,fc + 1}, and N^''+'^'^{t,^) e S^{-k - l,k + l} again fol- 
lows from Lemma [4.2| Furthermore, this choice of matrices implies Bk+i{t,S,) S 
5*{— fc — 1, fc + 2} and the recursion step is completed. 
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4.4.3. Zone constants and invertibility. If we consider only the part of the phase 
space defined by (1 + t)\£^\ > c for large c, the matrix-norm of I — Nk{t,(,) £ 
5*{— 1,1} is of size ©(c^^) and choosing c large enough implies invertibility of 



the matrix family Nk{t, ^) with iV^T (i, e 5*{0, 0}. Hence, the problem ( |4.27[ ) is 
equivalent to considering 

Dtl/W = {V{t,0 + Fk-i{t,0 + Rk{t,0)V^''^ (4.37) 

for y('=)(i,^) = N^:^{t,^)V^"\t,^) and with remainder term 

i?fc(i, = -Nk\t, ^Bkit, e 54-fc, fc + 1} (4.38) 

within the smaller hyperbolic zone Z\^yp{ck) = {{t,(,) : (1 + 1)\£^\ > Cfc}. 
We collect the main result in the following lemma. 

Lemma 4.3. Let k G N, k > 1. Then there exist matrices Nk{t,(,) E 5* {0,0}, 
diagonal matrices Fk^i(t,^) G 5* {0,1} and a remainder Rk G S^{—k,k + l}, such 
that 

{Dt-v{t,0-Ro{t,0)Nk{t,0 

- Nkit, (Dt - P(i, - Fk^t, - Rk{t, 0) (4.39) 

holds true within .^hyp(cfc) for sufficiently large zone constant depending on k. 
The matrices Ni^{t,£,) are uniformly invertible within Zi^yp{ck) with Nj^ it,£,) G 
540,0}. 

4.5. Solving the diagonalised system 



It remains to asymptotically solve (4.371 as (l+i)|CI — ^ 00. This will be done in two 
main steps. We always assume that we are indeed working inside the hyperbolic 
zone 2hyp(cfc) for sufficiently large constant c^. 

4.5.1. Treating the diagonal terms. The fundamental matrix of the diagonal part 



of (4.371, i.e., the matrix-valued solution £fe(t,s,^) of 

BtSkit, s, = {V{t, + Fk^iit, 0)£kit, s, 0, £k{s, s, = I e C^^^ (4.40) 

is just given by integration and taking exponentials 

£k{t,s,0^exp(ij^ {V{T,0+Fk-i{T,0)dTY (4.41) 

Here we essentially used that diagonal matrices commute with each other. By 
assumption, the entries of I?(t,^) are real and yield a unitary matrix after expo- 
nentiating. Similarly, i^(^^(t,^) G S^,{—j,j + 1} is integrable over the hyperbolic 
zone for all j > 1. This yields that the actual large time asymptotic behaviour of 
£k(t,s,£,) is encoded in the term F^^\t,^) and we can show that 



\£k{t,s,0\^exp(~l^ ImF(")(T,OdT 



(4.42) 



holds true for all individual diagonal entries of £k{t,s,£,) as two-sided estimate 
with constants depending on k but not on s, t and ^. From F''^\t,£_) G 5*{0, 1} 
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we can only conclude polynomial bounds on the right hand side, thus there exists 
exponents Ki and K2 such that for t > s the estimate 



1 



l + t 



Ki 



<\\£k{t,s,m 



< 



l + s 



K2 



(4.43) 



is valid within Zhyp(cfc). The exponents are independent of k and a similar estimate 
is true for t < s where the exponents are switched. 

4.5.2. Generalised energy conservation. We will speak of generalised energy con- 
servation for the system ( |4.15[ ) , if its fundamental solution (given as solution to) 

Bt£{t,s,O^A{t,0£(t,s,0, £-(5,5,0 = IeC'^^' (4-44) 

satisfies 

\\£{t,s,m<'^ (4.45) 

uniformly in it,£,),{s,^) £ Zhyp(c) for some c (regardless of their order). The 
generalised energy conservation property can be fully characterised by the term 
F^^H^iO — diagi?o(i,<^) constructed within the diagonalisation procedure. 



Theorem 4.4. The system (4.15) has the generalised energy conservation property 
if and only if 



sup 

(t,?):('5,«)e2hyp(c) 

holds true for some zone constant c. 



J s 



< 00 



(4.46) 



Sketch of proof. Since the matrices M{t, ^) G 5*{0, 0} and Nk{t, ^ G 5*{0, 0} have 
uniformly bounded inverses M~^{t, ^),N^^{t, £_) G 5*{0, 0}, the generalised energy 



conservation property of system (4.15) is equivalent to that of system (4.37). We 



exploit this fact in combination with the polynomial bounds (4.43) which are also 



valid (with the same exponents) for the fundamental matrix Skit, s,^) of (4.37). 

The main idea now is, that when choosing k large enough, the polynomial 
decay of the remainder term Rk (i, ^) compensates the increasing behaviour and 
we can use the Duhamel representation 

/oo 

£k{t,e,ORk{e,0£k{o,s,Ode (4.47) 

for k> Ki+ K2 with 

Zk{s,0=l + i £k{s,9,ORk{0,O£k{9,s,Ode (4.48) 



as integral equation relating £k{t, s, ^ to £k{t, s, C)- Indeed, 






1 



+ sj (i + e')fe+i \i + 



1 



de 



-k 



((l + ,s)|^|)-^ (4.49) 
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and similar for 



£k{t,e,ORkio,O£kie,s,Od0 



< 



Ki-k( 



i^r^(i + t)^i-'^-(i + s) 



\-K2 



(4.50) 



Therefore choosing k and s large implies that the difference £fc(t, s, ^) — £k{t^ s, ^) 
can be estimated in terms of (1 + t)^^{l + s)~^^~'^\£,\~'' (coming from Zk) and 
(1 + i)^i~'^(l + s)~^^ If I"'' (coming from the integral) and we can treat Skit, s, f ) 
as small perturbation of Skit, s, f ). 



The statement now follows from the observation that (4.46) is equivalent to 
a uniform bound on both Skit, s, f ) and its inverse. D 

It should be clear from the statement of the previous theorem, but we will in 



particular point out that the integral condition (4.46 1 is independent of the choice 



of the diagonaliser Af(t,^) of Ai(i,^) as long as it satisfies the symbol conditions 
Mit,C),M^^{t,^) e 5*{0,0}. Condition l\4A6l does depend on the lower order 



terms A(t,^) ~ Aiit,^) modulo 5*{— 1, 2} and on the time-evolution of the system 
of eigenspaces of v4i(t, ^). 

We will give an example. If we consider the damped wave equation 

utt - Au + 2b{t)ut = 0, (4.51) 

it can be rewritten as a system of first order in (|^|u, Dfu)^ with coefficient matrix 

M mt), 

Applying the diagonalisation scheme of Section |4.4| to this symbol yields terms 
Vit,i) = diag(|e|,-|CI) e 5,(1,0} and Fo(t,e) = ^^(i,^ = i&(i)I e 540,1}. 



Ait,i) = 



e 541,0}. 



(4.52) 



Hence, the generalised energy conservation property holds for (4.51 1 if and only if 
the integral 



6(r)dT 



(4.53) 



remains bounded. In case that 6(t) > 0, this just means h € L^(]R+) and solutions 
are even asymptotically free in this case. If the above integral is not uniformly 
bounded, we obtain the two-sided estimate 

\Skit,s,0\~e^v(j^biT)ATY (4.54) 

This is consistent with the observations of Section 13.1.21 for the scale invariant 



model case. The equation (4.51) was studied by the second author in [65] under 
similar assumptions. 



4.5.3. Perturbation series arguments. For the following we assume that (4.46) is 
valid such that in particular 

||£4t,s,0ll«l (4.55) 

for all fc > 1 uniformly on (i,f), (s,0 ^ 2^hyp(cfe). Then we can follow a simpler 
argument to construct Sit,s,S,) explicitly. 
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We make the ansatz 

f/c(i,s,OS/c(i,s,0 (4.56) 

for the fundamental solution of the diagonalised system Dt — P — F^-i — Rk with 
a so far unkown matrix function Qk{t, s,^). A short calculation shows that the 
function Qk{t,s,£,) solves 

BtQk{t,s,0=nkit,s,OQk{t,s,0, Qk{s,s,0=l, (4.57) 

with 

TZkit, s, = £k{s, t, ORk{t, 0£k{t, s, 0- (4.58) 



We can represent the solution to (4.57) in terms of the Peano-Baker series 

Qk{t,s,o = i + J2i' nk{ti,s,o nk{t2,s,o 

rti-1 

••• / nk{ti,s,(,)dU---<lt2dh. (4.59) 

J s 

If fc > 1, the symbol estimate Rk{t,S,) G S^{—k,k + l} allows to show the bounded- 
ness of Qfe(i, s,^). Indeed, using \\TZk{t,s,(_)\\ w ||i?fe(i,f)|| uniform in s we obtain 

\\Qk{t,s,0\\<expU \\nk{T,s,0\\dT)<exp{Cc-'') (4.60) 

uniformly in Zhyp(c). 

We will derive one more consequence from the series representation. The 
uniform integrability of TZk{t, s,^) over 2hyp(cfc) implies that the limit 

lim Qfc(i,s,0 = Qfc(c^,s,0 (4.61) 

exists locally uniform in (5, ^) G 2^hyp(cA:)- This is most easily seen by checking the 
Cauchy criterion or reducing it to the estimate 

/oo 
\\Rk{T,m\\Qk{r,s,0\\dT 

/OO 
II i?,(r, oil dr. (4.62) 

The matrices Qk{t,s,£,) as well as Qfc(oo,s,0 are invertible as the uniform lower 
bound on their determinant 

det Qfc(i, s,0 — exp I i / trace7?.fc(r, s,Odr J 

> exp (~d I ||i?fc(T,OI|dT j > C> (4.63) 

shows. 

In the sequel we will also need estimates for ^derivatives of Qk{t, s, £,), which 
can be obtained by differentiating the above representations provided k is sufB- 
ciently large. There are some subtleties involved. If we consider ^derivatives of 
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£fc(i, s,^) then logarithmic terms might appear from integrals of ^-derivatives of 
F<-°\t,0 e S40,l}, i.e., we get 

l!Dfffc(i,ie'OII<(l + 0'"'(log(° + 0)'"' (4-64) 

for t^ solving (l + ij)!^] ~ Ck (or being equal to if this solution becomes negative). 



This implies that derivatives of 7?,fe(i, s,^) defined in (4.581 satisfy weaker estimates 
than the symbol estimates of Rk{t,^). Based on the definition of the hyperbolic 
zone and the above estimate we obtain 

llD^7^fe(t,i5,OII<ler'-|"l(l + t)"'(log(e + i))'"' (4.65) 

for \a\ < (fc— 1)/2 and uniform in Zhyp(cfc). Differentiating the series representation 



(4.59) term by term and using this estimate yields 

\\^tQk{t,^,m < i^ri"i(iog(c+^))'"' (4-66) 

uniform in 2^hyp(cfc) and with constants depending on fc. 

4.6. Examples and resulting representations of solutions 

We will discuss some examples in more detail in order to show how to apply 
the previously developed theory. We will restrict our attention to homogeneous 
problems as in this case the estimates within Zp^ are easily obtained. In fact, 
assuming that A{t, ^) is positively homogeneous of order one in ^ yields \\A{t, ^)|| < 
C|^| uniform in t and, therefore, also directly the estimate 

||£(t, 0,011 < exp (J \\A{r,m<ir] < exp(Ct|^|) < 1, {t,0 e Zpd, (4.67) 



for the fundamental solution defined by (4.44). Furthermore, differentiating the 



differential equation (4.44) implies that derivatives satisfy 

l|Df£:(tc,0,OI!<a|frl"l. (4.68) 

In the general case, the treatment for small frequencies needs more care and is 
often reduced to solving Volterra equations in suitable weighted spaces, see [65] . 
[U or 0. 

4.6.1. Symmetric hyperbolic systems. First, we consider symmetric hyperbolic 
differential systems of the form (|4.1|) , 



BtU - ^ A,(t)D,^.{/, (7(0, •) - Uo, (4.69) 

for self-adjoint matrices Aj{t) E T{0}. By partial Fourier transform this is related 
to DtU = A{t,£,)U with 

n 

A{t,0 ^ Y. Mm ^S4l,0}. (4.70) 
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We assume that the roots of ( |4.4[ ) satisfy 

|A.(i,e)-A,(t,e)i>ciei 

uniformly in t and i =/= j for some constant C > 0. Then we can follow the scheme 



of Sections 4.4 and 4.5 to explicitly determine leading terms of the representation 
of solutions. 



As (4.701 defines a self-adjoint matrix, we choose the diagonaliser M{t,^) 
as unitary matrix depending smoothly on t and ^. This gives M{tX) G 5* {0,0}. 
Furthermore, denoting the columns of the matrix M(t, f ) as Vj{t, ^), j = 1, . . . , d, 
we can express the first matrices constructed in the diagonalisation scheme as 
I?(i,0 = diag(Ai(t,0,---,Ad(t,0) and 



(i^W(i,0),, = -i9*^,(t,0 • ^,(i,e) - -\dt\\v,{t,0\? - 0. 



(4.71) 



The next diagonal terms starting with i^'^^)(i,^) are in general non-zero such that 
-Ffc_i(i,^) e 5*{— 1,2} in place of the weaker result from Lemma 4.3 In conse- 



quence, the estimates (4.66) are valid without the logarithmic terms. 



Combining the representation (4.561 with the representation ( 4.41[ ) and the 
estimates (4.66) and using the symbol properties of the diagonaliser M[t,^) and 
Nk{t,£,) together with estimate (4.68) we obtain structural information about the 



representation of solutions. Note first that for any t > t^, i.e., within the hyperbolic 
zone Zhyp(cfe), the representation 

7Vfc(i4,C)Af(tc,O^(i?,0,O (4.72) 

holds true. Regrouping the expressions a bit, we obtain the following theorem. Note 
that the precise zone constants depend on the number of derivatives of Bj(t,^) we 
have to estimate. 



Theorem 4.5. Any solution to the system (4.69) is representable in the form 

uit,0 



d 

E 



e'*''^(*'«)i?,(i,Of/o(e)+fpd(t,e)C/o(0, 



(4.73) 



where £pd(t,^) is uniformly bounded and supported within the zone Zpd(2cfe), the 
phase function is defined in terms of the homogeneous characteristic roots Xj{t,^) 
via 

Mt,0^]J X,ir,Odr (4.74) 

and the amplitudes Bj{t,^) are supported within Zhyp(cfc) and satisfy the symbolic 
estimates 

for all \a\ < (fc - l)/2. 



(4.75) 
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We omit an explicit proof as it is just a combination of the above men- 
tioned estimates and representations. We only point out that one uses a variant 
of Xhypit,£,) to cut off the zone Zpd{ck) and to decompose the representation 
smoothly into parts. 

4.6.2. Second order equations. We consider first homogeneous second order equa- 
tions of the form 

n 

dtu - ^ a,j{t)dx,dx^u = 0, u(0, •) = uo, wt(0, •) = "i, (4-76) 

for aij{t) g 7~{0}. The assumption of uniform strict hyperbolicity is equivalent to 
assuming that 

n 
J2 <^^Am^^, > C\e (4.77) 

for some constant C > 0. Let a(i,0 denote the (positive) square root of the 
left-hand side of this inequality. Then a(i,^) G 5*{1,0} and if we denote U = 



{a{t,£^)u,Dtu) , then equation (4.76) is reduced to the first order system 



_ / BMW , .A_ 

D,C/ = A{t,OU = f/) "^''^^ [/, [/(O, •) = t/o. (4.78) 

\a(t,Q J 

The above equation and this system were treated by Reissig |45| and in a similar 
setting Reissig- Yagdjian [411 . 

Clearly A(i,^) e 5*{1,0}. Furthermore, a diagonaliser of the homogeneous 
principal part is given by 

M-{\ ;), M-'.i{\ ;), ,.™, 

so that V{t,^) = diag(a(t,C),-a(t,0) and F^°\t,^) = ^^^I. The matrix 
£o{t,s,^) can be calculated explicitly as 



£o{t,s,0^exp(ij^ (l?(T,0 + F^°'(T,0)dT 



ya(^^. 



^ diag f exp f ±i / a(T,^)dT 



(4.80) 



Va(s,0 
As can be seen clearly from this representation, the exponential of the primitive of 



F^^H^jO is a symbol from 5* {0,0} and the logarithmic terms appearing in (4.66) 



do not occur. Furthermore, estimate (4.67) applies for the fundamental solution 



associated to {u,Dtu)^ and by equivalence of norms also to the one of (4.78). 



Hence, we obtain a representation of solutions reminiscent to that of Theorem|4.5| 
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Theorem 4.6. Any solution to (4.761 can be written as 






"(^'0 = E Ee^"'''*'^^^^-.±.hyp(i,0 + fe..pd(t,0 u,{0 (4.81) 



with 6j.pd(i,C) uniformly bounded and supported within Zpd(2cfe), the phase func- 
tion given in terms ofa{t,^) as 

^(^,0-7 / «(T,Odr (4.82) 

^ Jo 

and amplitudes fcj.±,hyp(t,0 supported within Zhyp(cfc) and satisfy 

|D?^.,±,hyp(t,OI<a|Crl"l-^' (4.83) 

/or |a|<(fc-l)/2. 

4.7. Dispersive estimates 

In order to derive dispersive estimates from the results obtained so far, we have 
to provide a suitable parameter dependent version of the appearing oscillatory 
integrals. Such parameter dependent versions were developed in the setting of 
equations of constant coefficients with lower order terms in }49| based on the 
multi-dimensional van der Corput lemma in [^ . For the presentation below we 
foUow WD and ^. 

4.7.1. Contact indices for families of surfaces. In order to formulate estimates in 
a uniform way over a family of surfaces, we need some more notation generalising 
the treatment of Section |1.3| Let again S be any closed smooth hypersurface in 
M" and let p G S be any point. We translate and rotate the hypersurface in such 
a way that it can be parameterised as 

{{y,h{y)):yen} (4.84) 

in a neighbourhood of the point p for a suitable open set il C M"~^. We start 
by considering the convex situation. 2-planes H containing the normal NpTi are 
determined by directions w S S"~^ in J7 C M.^~^ and thus the contact orders 
7(E;p, iJ) correspond to the vanishing order of the function p n- h{puj) at p = 0. 
We need to measure this vanishing in a more quantitative way and define for some 
given parameter 7 S N>2 

7 
^i^,r,p) = inf E |^^^(P^)lp=o| ■ (4.85) 



From the definition of the contact index 7(S) in ( 1.18 1, it immediately follows that 
>r(I], 7(I]);p) > for all p G E. Because the function p i— > >f(I],7;p) is continuous 
and E closed by assumption, the minimum 

k{J:, 7) = min >c{j:, 7; p) (4.86) 

pes 

exists and is also strictly positive for 7 = 7(E). 
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If we find a uniform lower bound on >c(T,x, 7) for a family of convex surfaces 
Sa depending on parameters A G T for some number 7 > 2 then we say that 
this family satisfies a uniformity condition. We further define the uniform contact 
index 

7unif ({Sa : a e T}) - min{7 e N>2 : inf k{^x, 7) > 0}. (4.87) 

AG 1 

For real parameters t e M-)_ (which can be thought of as time for now) it makes 
sense to consider a weaker notion of asymptotic contact index and we define 

7as({St : t > to}) - min{7 G N>2 : liminf ^(S^, 7) > 0}. (4.88) 

In the particular case that T,\ is a family of spheres of different radii r^, the 
constant >r(I]A,2) « r^^ is a measure of the curvature of the sphere and we 
have 7unif ({Sa}) = 2 if and only if rx stays bounded. If on the other hand r\ is 
unbounded, the uniformity condition is not satisfied and no uniform contact order 
exists. 

For non-convex surfaces we need to define a similar uniform non-convex con- 
tact index. We only explain the differences to the above formulae. First, we replace 



(4.85) by 



7 



M^,T,P) = sup Y, \d'ph{puj)\p^o\ ■ (4.89) 

l'^l = lj=2 

and then set kqCE,^) = imiip^-s kq{Ti,j;p). Again >fo(S,7o(S)) > and we use 



the analogues to (4.87) and (4.88) to define uniform and asymptotic contact orders. 
If a family Sf satisfies the uniformity condition, then the constants in the 
estimates of Lemmata 1.3 and|1.5|are uniform over the family of surfaces. 



4.7.2. Estimates for t-dependent Fourier integrals. Now we are in a position to 



discuss Fourier integrals appearing in the representations (4.73) and (4.81), or, 
more generally, representations obtained by the construction of Section~ [45l We 
omit the possibly occurring logarithmic terms here; they can be included easily 
with a slight change in the decay rates. 

We consider a i-dependent family of operators 

Tf.uo^ [ e'(--«+*''(*^«))a(t,0"o(OdC (4.90) 



for a real- valued homogeneous phase function 'd{t,£,) G C°°(M+ x (M" \ {0})) 
satisfying 'd{t,pS,) — p^(i,^) for p > and 

C-i|e|<^(i,0<qCI, |D5"i?(t,ai<a|Cr"l"l, (4.91) 

for some constants C > 0, Ca and all t > tg, and an amplitude a(t,5) supported 
within Zhyp(c) and satisfying the symbol estimates 

|D^a(i,OI<a|erl"l (4.92) 
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for a certain (finite) number of derivatives. Associated to the phase function we 
consider the family of slowness or Fresnel surfaces 

S4 = UeM":79(t,e) = l}. (4.93) 

Then the following theorem can be obtained by the method of [5T| . 

Theorem 4.7. Assume the family Sf is convex for t > to, satisfies the uniformity 
condition and that^ — 7as({St}) is its asymptotic contact index. Then the estimate 

\\Tt\\Bi,^L^<Ct-^ (4.94) 

holds true for r = n — ^^^^ . 

We remark that only finitely many derivatives of the phase (meaning also 
finite smoothness of St) and of the symbol a{t,^) are needed to prove this state- 
ment. 



Proof. Using the definition of the Besov spaces from (1.23) together with Holder 
inequality we see that it suffices to prove 

\\Ttoq^{2-^\B\)uo\\L-<2^^''-'^h-^\\uo\\L^, (4.95) 

i.e., corresponding bounds for the operator with amplitude a(t,^)(/)(2^-'|^|) taking 
into account the desired Besov regularity. The uniform bound (|4.91) on the phase 



allows us to find a function V G Cff (M+) satisfying 0(|^|)V'(i9(i,O) = '/'(ICD- We 
denote by 

Ij{t,x)^ /e'("-«+*^(*'«»a(t,0V'(2--''t?(i,e))dC (4.96) 



the (smooth!) convolution kernel of an operator related to (4.95 ) and it is sufficient 
to prove the uniform bound 

sup \Ij{t,x)\ < 2^(""'^'r'^, j > 1, (4.97) 

a:GE" 

in combination with a related low-frequency estimate. The bound in Ij{t, x) can be 
achieved by the stationary phase method in combination with the multi-dimensional 
van der Corput lemma due to |48j . 

Stationary points of the phase are solutions to x + tV^^(i,^) — 0. We use 
a cut-off function x G C§°{M.") with x(a;) = 1 for small |a;| and decompose the 



integral (4.96) into two terms 



I^'\t,x) = j e'(-«+*''(*^«))x(:^/i + V5z9(i,C))a(i,e)V^(2"^,?(t,e))de, (4.98) 

I^^\t,x) = /e'(-«+*''(*'«»(l-x(a;A + V^z?(t,0))a(i,0VX2-^"^(i,0)de (4.99) 



Based on |a;-|-iV^^(t, ^)| > t, the second integral can be treated by integration by 
parts giving 

1/^(^,01 < CArt-^2^("-^) (4.100) 
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uniform in j and t and for any number Nr\ For the first integral we use the 
strueture of the level sets Sf and restrict consideration to large values of t. We 
localise the integral into narrow cones, by translation and rotation we can assume 
that ^ is within a sufficiently small conical neighbourhood of (0, ... , 0, 1)^ and St 
can be parameterised as {{y,ht{y)) : y £ U} with h vanishing to at least second 
order in 0. The function V/it : U — >■ Vht{U) C M"~^ is a homeomorphism and ht 
is concave as Ej was assumed to be convex. Furthermore, we have 

|D«/it(y)|<a, t>io, (4.101) 

uniform on [/ if [/ is small enough. Indeed, as '&{t,y,ht{y)) = 1 we obtain by 
differentiation Vj,i9 + d^^dVht — and by Euler's identity d^^'d{t, e„) = d{t, en). 
Therefore, the bound |Vyi?(i, y, ft.)| ^ 1 implies the bound on \/ht{y) uniform in t. 
Higher order derivatives follow in analogy. 

Associated to Sj we have the Gauss map 

and for given x ~ (x',x„) near — iV^i?(t, e„) we define Zt G C/ by G{zt,ht{zt)) — 
—x/\x\. Making the change of variables ^ = {ry,rht{y)) with i9(t,^) — r Ki 2^, we 
get that the localised part of /■ {t,x) equals 

M.'-y+.„h,{y)+t)^(^^^ ry, rht{y))i:{2-^r)x{t, x, y)K{t, r, y)dydr 

= 2^"/ / e'2'''(^'2'+^"'''(!')+*)aj(i,r,y)x(t,2:,y)d2/dr (4.103) 
Jo Ju 

with x{t,x,y) ^ x{x/t + V ^d{t,y,ht{y)), K{t,r,y) = \dS,/d{r,y)\ the Jacobi deter- 
minant of the i-dependent change of variables, and 

&,it,r,y) = ii^{r)a{t,Vry,Vrht{y))K{Vr,y) (4.104) 

satisfying uniform bounds of the form 

\D'^~a{t,r,y)\<C^. (4.105) 



Here we used the symbol estimates of a in combination with (4.101) and also 
that the latter imply uniform bounds on y-derivatives of the Jacobian. Note, that 
r-derivatives of the Jacobian have symbolic behaviour by homogeneity. The sta- 
tionary phase estimate i48j Theorem 2.1] is applied to the integral over t/, which 
gives 



^^'-i-'-y+-^^^iv)+^)~a^{t, r, y)x{t, x, y)dy 



< C\Vxn\~^^ (4.106) 

with a constant C uniform in the remaining variables. Due to the localisation we 



have that |a;| w |a;„| w t and therefore we obtain the desired bound (4.95) 



^Of course, A^ depends on the number of derivatives we can estimate by ( |4.92[ l. The number that 
we need is [(n — l)/7] . 
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The remaining low-frequency estimate follows in an analogous manner, omit- 
ting the dyadic decomposition and using the relations t\^\ > 1 combined with 
the symbolic estimates and the idea that uniformly bounded Fourier multipliers 
supported within Zpd always lead to t~" decay rates (which are much faster). D 

If the convexity assumption is dropped, decay rates get much worse but on 
the other hand we also need less regularity in the estimate. We only formulate the 
theorem. 

Theorem 4.8. Assume the family Et satisfies the uniformity condition and that 
7o — 7o,as({St}) is its non-convex asymptotic contact index. Then the estimate 

||rjB._^^Loo <cr^ (4.107) 

holds true for r = n . 

4.7.3. Extensions to fully variable setting. We conclude Section |4] with some more 
general estimates. They apply to global Fourier integral operators of a very par- 
ticular structure at infinity and are essentially based on the observation that the 
variables t and x appear parametric in dispersive estimates. They also fit to the 
earlier observation that i^ is a natural co-variable when describing solutions to 
free wave equations and equations with weak dissipation. 
We consider the operator 

Tf.uo^ f e'(^-«+*''(*'^'«»a(t,a;,e)uo(Ode (4.108) 



for a homogeneous real phase function ??(t, x, £,) € C°°(M+ x M" x (M" \ {0})) with 
'd{t, X, p^) = p'd{t, X, ^) for p > and 

C-'\^\<^t,x,0<C\^\, |D^",9(i,x,OI<ier"'"', (4.109) 

and the amplitude a{t, x, S,) supported in {(i, x, ^) : (l+i)|C| > c} for some constant 
c and subject to 

|D^"a(i,x,OI<C7„|C|-l"l. (4.110) 

A major difference to the previous considerations is that we now have a family of 
surfaces parameterised by t and x 

St,. = UeM":i9(t,a;,0 = l} (4.111) 

and we have to use a contact index defined asymptotic with respect to t and 
uniform in x. Let therefore 7 = 7as,unif ({^^t.a;}) with 

7as,unif({St,a.}) = min{7 G N>2 : liminf inf >f(St,.,7) > 0} (4.112) 

— t— i-oo xGR" 



and we assume that all the surfaces are convex. Then the proof of Theorem |4.7| 
carries over and gives 

\\TM\L^<t-^\\u^\\Br^^ (4.113) 
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Tt-l 

7 

the weaker decay rate f^/'^o for the non-convex asymptotic uniform contact index 



for r > n — - — . A similar replacement works for the non-convex situation with 



defined in analogy to (4.112) 



4.8. An alternative low-regularity approach: asymptotic integration 

We will conclude Section |4] with a short review and reformulation of recent results 
due to Matsuyama- Ruzhansky 30J . which allow for very low regularity of the 
coefficients. The approach originates in asymptotic integration arguments due to 
Wintner [61j. The price to pay is that results are not optimal with respect to 
regularity and they do also not reach the critical threshold for lower order terms. 
We will formulate them for hyperbolic systems as it is in the context of the 
present exposition. The case of scalar equations has been treated in |30j while 
the case of systems can be found in i31j . To be precise, we consider the Cauchy 
problem 

DtU ^ A{t,r>)U, f/(0,-) = t/o, (4.114) 

for initial data Uq e ^^(M";^*) and with symbol split as A{t,^) = Ai(i, C) + 
Aq (i, ^) into a homogeneous strictly hyperbolic part 

A,it,p0^pA,it,0, Vp>0, (4.115) 

having real eigenvalues specv4i(t, ^) = {Ai(t, ^), . . . ,\d{t,^)} C M uniformly sepa- 

inf '^'(^'^)-^^-(^'^)' > (4.116) 

for all i ^ j. Difference to the previous considerations is that we do not assume 
smoothness with respect to t, instead we assume that the symbol is of bounded 
variation and satisfies 

\^\~'dtAii;0 e L\R+;&'"'), Aoi;0 e L\R+;C'''"') (4.117) 

uniform in ^ 7^ 0. 

These assumptions have some almost immediate consequences. 

Lemma 4.9. 1. \£_\-'^dtXj{t,£,) £ L^{R+) uniform m£_^0. 
2. There exists a matrix N{t,^) diagonalising Ai{t,S,) 

^i(i,OA^(t,0 = ^(i,Odiag(Ai(i,0,...,A<i(t,0), C 7^ 0, (4.118) 

with \\N{t,0\\,\\N^^it,0\\ uniformly bounded and dtN{t,^) £ L\R+) uni- 
form in ^ ^ 0. 

Sketch of proof . The second statement is analogous to [33l Prop. 6.4] and also to 
Lemma 4.2 3, while the first statement follows from differentiating (|4.118) with 



respect to t. D 

We can use N{t, ^) to transform the system into a diagonal dominated form. 
Considering V{t,£^) — N^^{t,£^)U{t,^) we obtain the equivalent system 

BtV={V{t,0 + Rit,0)V, (4.119) 
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with P(t,e) = diag {Xi{t,0, ■■■, Xd{t,0) and 

Rit,0 = N-Ht,OMt,ONit,0 + {BtN-\t,0)N{t,0. (4.120) 

By assumption on Ao{t,£^) and construction of N{t,^) it follows that i?(-,0 ^ 
L^(]R+) uniform in ^ and we can solve the resulting system by an argument similar 
to Section |4.5| The only difference is that we obtain just uniform results with 
respect to S, and have no information about the behaviour of derivatives with 
respect to ^. 



Theorem 4.10. Solutions to (4.114) can be written as 



U{t,0 = Ee"'''^*'^^^^(^'^)^o(^) (4.121) 

with phase functions given in terms of the characteristic roots Xj{t, ^), 

^jit,0 = \J ^,iO,OdO, (4.122) 

and amplitudes Bj{t,S,) G L°°(M+ x M") satisfying 

B,{t,0 = a,(0 + £,(t,0, lim e,{t,0 = (4.123) 

locally uniform in ^ ^ 0. 

If one has additional information on the behaviour of ^-derivatives of the 
symbols ^^(i,^), k = 1,2, one can prove better estimates and in particular also 
obtain symbolic type behaviour of the amplitudes Bj{t,^). We give only one cor- 
responding statement. If Ao{t,^) — and Ai{t,^) satisfies 

\(\\"\''D^dtAi;OeL\R+) (4.124) 

uniformly in ^ 7^ 0, then one can show that 

||D^i?,(i,OII<a, Cy^O, (4.125) 



and dispersive estimates based on Theorems 4.7 and 4.8 the (explicitly known) 
finite number of derivatives needed there can be obtained by requiring a higher 
regularity of initial data. We refer to [3ni [U] for further details. 



5. Effective lower order perturbations 

If lower order terms are too large to be controlled, it becomes important to investi- 
gate the behaviour of solutions for bounded frequencies. We will restrict ourselves 
to situations where an asymptotic construction for ^ — ^ becomes important and 
provide some essential estimates for this. 
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5.1. The diffusion phenomenon 

The classical diffusion phenomenon gives an asymptotic equivalence of damped 
wave equations and the heat equation. It was first observed for porous media type 
equations and its formulation for the damped wave equation is due to Nishihara 
[55] . P^ and independently Han-Milani [32] in various formulations. The esti- 
mates of Nishihara were extended to arbitrary dimensions by Narazaki |34| and 
provide U'-L'^ type estimates for differences of solutions to damped wave and heat 
equation. 

We will follow a different line of thought here and provide estimates in the 
flavour of Radu-Todorova-Yordanov [38]. They are only energy type estimates 
and therefore easier to obtain. On the other hand, they are flexible enough to be 
formulated for abstract Cauchy problems in Hilbert spaces based on the continuous 
spectral calculus for self-adjoint operators. 

Let in the following u{t,x) be the solution to the Cauchy problem for the 
damped wave equation 

Utt-Au + Mt=0, ^(0, O^wo, ut{0,-)=ui, (5.1) 

and similarly v(t,x) be the solution to the Cauchy problem for the heat equation 

Vt = Av, v{0, ■) ^vq = uo + ui (5.2) 

with related data. Then their solutions are related due to the following theorem. 
Note that in general 

Mt, ■)\\l2 < C{\\uo\\l2 + \\ui\\h-i), Ht, ■)\\l- < C\\vo\\l2 (5.3) 
are the best possible (operator-norm) estimates for solutions of both problems. 



This can be seen similar to the discussion in Section 1.2 using explicit represen- 
tations of Fourier multipliers. However, forming the difference of the solutions 
improves estimates by one order. We denote by e* the heat semigroup. 



Theorem 5.1. Let u(t,x) and v{t,x) be solutions to (5.11 and (5.2|, respectively. 
Then 



1. the difference of the solutions satisfies 

\Ht, ■) - v{t, ■)U. < t-i (||e*^/2^olU^ + \\c'^/^u,U.) 



+ e-*/i«(||uo|U. + ||i.i||^-i) (5.4) 

for all t > 1; 
2. moreover, for all k (z N and a G Nq the higher order estimate 

||D?DS(«(i, .) - v{t, ■))U. < i-i-'=-l"l/2 {\\e'^/'u,\\,2 + |le*^/V|U^) 

+ e-*/iS(||Mo||H<=+i«i +||ui||fffc+ic|-i) (5.5) 
holds true for t > 1. 
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Proof. As the theorem extends the results of [3.8^ we include a sketch of the proof. 
Fourier transform relates the Cauchy problems to initial value problems for or- 
dinary differetial equations which can be solved explicitly. The solution of the 
parabolic problem v{t, £,) is given by 

v{t,0^e-'^i^\uo{0 + MO), (5.6) 

while the solution to the damped wave equation can be written as 

u{t,0 = e-*/^(cos(i^^^^T7i)So + '-^^^^^^f^^MO) (5.7) 

VKr -1/4 

for 1^1 > 2 and similarly 



v/V4-iep 



(5.8) 



for 1^1 < 4 As we are concerned with polynomial estimates, we can neglect all ex- 
ponentially decaying terms (and collect them in the e~*' ^^ estimate later on). This 
happens for v{t,^) if |^| > 1/4, for u{t,^) for |^| > 1/2 and for the exponentially 
decaying terms defining cosh and sinh. What we are left with are the terms 

v{t,x): e-*l«l'(uo(0+Si(e)) (5.9) 

u{t,x): c-*/2WV4-|«P (Co(Ouo(0 + CiiOuiiO) (5.10) 

locally near ^ = and with radially symmetric smooth functions Co{S,) and Ci(^) 
satisfying Cq{0) = C'i(O) = 1. Our first observation concerns the exponent in 



(5.10). It satisfies 



v/l/4-|cp = l/2-|el'-Vfefclcp^ 6a. >o, (5.11) 



fe=2 



so that the main terms of the difference u{t, ^) — v{t, ^) are a sum of terms of the 
form 

e-*l«P (^^-tb.\i\'+toi\if) _ i) (S„(^) + u^iO), (5.12) 

e-*l?l^+*o(|?r)((C'o(0 - l)uo(0 + (Ci(0 - l)MO)- (5.13) 

The first statement is proven, if we can bound these two terms by the ones on the 
right hand side; i.e., we are looking for uniform bounds for the multipliers 

te-*i?iV2 (^e-b.m'+tom'^) _ i) (5.14) 

and 

te-*l«l'/2+*o(l«l')(C,(e)-l), J =0,1. (5.15) 

Using the elementary estimates |e^* — 1| < s and se~* < 1 for s > we conclude 
that the first multiplier is bounded by t^|^|'*exp(— i|^p/2) < 1, while the second 
multiplier can be estimated by j^pte"*'^' < 1 and the desired statement follows. 
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For the second statement we observe two things. First, when estimating x- 
derivatives this just gives additional factors of ^ for smah frequencies and the above 
estimates improve by corresponding t~^''^ factors for each derivative. Similarly, 
when considering t-derivatives the representations get one additional factor of |^p 
for small ^ and each i-derivative improves the estimate by a factor of t~^. D 

5.2. Diagonalisation for small frequencies 

Versions of the diffusion phenomenon can also be obtained in a variable coefficient 
setting. We will restrict our consideration here mainly to i-dependent hyperbolic 
systems, for the sake of brevity even differential hyperbolic systems of the form 

n 

Dt[/ = ^Afc(t)D,,C/ + iS(OC/, U{t,-) = U^. (5.16) 

k = l 

Main difference to the considerations in Section [4] is that we assume now that B{t) 
is not of lower order in the T-hierarchy, i.e., we assume that 

Ak{t),B{t)(.T{Q). (5.17) 

We will make three main assumptions here; 

(Bl): the matrices Ak{t) are self-adjoint and B{t) > 0; 

(B2): the matrix B{t) has d eigenvalues = (5i(t) < (52(t) < • • • < Sd{t) satisfy- 
ing 

liminf |(5,(t) - 6Jt)\ > 0, i^ j; (5.18) 

(B3): the matrices A{t,^) = X^Li Ak{t)^k and B{t) satisfy for all u G C^ 

^ d-l 

-\\vr<j2e,\\B{t){A{t,oyvr<c\\vr, t>^, (5.19) 

for any choice of numbers eg, ... , e^-i > and suitable constants c and t^ 
depending on them. 

Assumption (Bl) guarantees that the system is symmetric hyperbolic and (par- 
tially) dissipative. Therefore, the energy estimate 

||C/(t,-)l|L^<l|C/o||L^ (5.20) 

is valid. By Assumption (B2) we know that one mode is not dissipated, while 
Assumption (B3) will be used to show that the high frequency parts of solutions 
are still exponentially decaying. Inspired by Beauchard-Zuazua [3 , we will refer 
to (B3) as uniform Kalman rank condition. If Ak and B are independent of t it 
just means that 

rank {B\A{S,)B\ ■ ■ ■ \A{Ci'^-^B) = d, (5.21) 

which is the classical Kalman rank condition arising in the control theory of or- 
dinary differential systems. Under certain natural assumptions, this is equivalent 
to the algebraic condition of Kawashima-Shizuta [H], but the latter are more 
complicated to rewrite uniformly depending on parameters. 

Our aim is to understand the large-time behaviour of small frequencies here. 
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Figure 2. Zones for the small- frequency diagonalisation 



For small ^ we use a diagonalisation procedure similar to the one from Sec- 
tion 4.4 It works within the elliptic zone 2on(c) = {(^,0 • ICI < c, i > c^^} 
collecting small frequencies for large times. All considerations will be done for 
t > to with to sufficiently large. This ensures uniform diagonalisability of B{t). 

5.2.1. Initial step. We denote by M{t) a diagonaliser oi B{t) satisfying the uniform 
bounds M{t),Nr^{t) e T{0}. Denoting V^°\t,0 = M-\t)U{t,^), we obtain the 
equivalent system 



Dty(") = iv{t) + Y,M-\t)Ak{t)M{t)a + {DtM-\t))M{t)\V^°^ (5.22) 
with T>{t) = idiag((5i(i), . . . , Sd{t)). We denote the two non-diagonal terms by 

n 

Ri{t,0 = ^M-i(t)Afe(t)Af(t)efc -I- {DtM-\t))M{t), (5.23) 



fc=i 



remarking in particular that the first term is linear in ^ with coefficients in T{0}, 
while the second one is independent of ^, but in the better class T{1}. Both are 
small compared to the difference 5i(t) — 6j{t) provided |^| is small and t is large. 
In the sequel we will use the notation 



r{m} = \ p{t,C) = J2 P'-WC ■■ Po'it) e n^ - \a\} 

\a.\<m 



(5.24) 
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for polynomials with coefficients in the T-classes. By construction, Ri{t, ^) e "Pjl}- 
5.2.2. The diagonalisation hierarchy. We start with a system 

Dtl/(") = (P(t) + i?i(t,e))y(°) (5.25) 

with Pi{t,£_) £ "Pjl}- Befor e set ting up the complete hierarchy, we will discuss its 



first step. Similar to Section 4.4 we construct a matrix A^'^-'^^ (t, ^) e Pll} such that 

(D,-P(t)~i?i(t,0)(I + iV(i)(t,0) 

~{l + M^\t,0){Dt^V{t,0-F,{t,C))eV{2} (5.26) 

holds true for some diagonal matrix Fi(t,^) G 7^{1}. Collecting all terms not 
belonging to the right class yields again conditions for the matrices N'^^'>{t,^) and 
Fi{t,^). Indeed, 

[2?(t),iV(i)(i,0]--i?i(i,0+^i(i,0 (5.27) 

must be satisfied and, therefore, we have 

F,{t,0 = diagi?i(t,0, [N('\t,0)^^ = l^t)_sl^ty (5-28) 

while we may again choose diagonal entries to be {N^-^'>{t,£^))i^i = 0. As desired, 
this implies iV(i)(i,0,^i(i,C) & ^{1}- 

Recursively, we will construct matrices A^*^''-' (i, ^) G V{k} and F^''^{t,Q G 
V{k} diagonal, such that for 

K K 

iVK(i,e) = I + E^'"(^'^)' FK{t,(,) = Y,F('^\t,0, (5.29) 

k=l k=l 

the estimate 

BK{t,0 = {^t-V{t) - R^{t,S.))NK{t,0 

- NK{t,0{^t-V{t) - FK{t,0) (^ViK + 1} (5.30) 

is valid. 

We just did this for X = 1, it remains to do the recursion k ^^ k + 1. Assume 
Bk{t, ^) G V{k + 1}. The requirement to be met is that 

Bk+i{t,0-Bk{t,0 = -[V{t),N<'^+^\t,S,)]+F^''+^\t,0 n,odV{k + 2} (5.31) 

for diagonal F''^+^\t,£,), which yields F''^+'^\t,i) = -diagBkit,^) together with 

Again, the diagonal terms can be set to zero, iN'^''^^^t,S,))ij — 0. It is evident 
that the construction implies F'^''+^\t, ^), N^''+'^\t, ^) G Vik + 1} together with 
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The matrices Nkit,^) e T'jO} are invertible with inverse N^^{t,^) e V{0} 
if we restrict our consideration to a sufficiently small elliptic zone Zcn(cfe). The 
result of the above consideration can be summarised in the following lemma. 

Lemma 5.2. Assume (Bl) and (B2) and let k E N, k > 1. Then there ex- 
ists a constant Cfe and matrices iVfe(i,^) S ^{0}, diagonal Fk(t,S,) G ^{1} o,iT-d 
Rk+i{t,0 e V{k + 1} such that 

(Dt-p(f)-i?i(t,e))iVfe(t,e) = ^fc(t,o(Dt-2?(i)-i^fe(t,e)-i?fc+i(t,e)) (5.33) 

holds true within Zou(cfe). The matrix Nk{t,S,) is uniformly invertible within this 
zone. 

It is worth having a closer look at the upper-left corner entry of the matrix 
Fk{t,£,) and consequences for them based on Assumption (B3). The above lemma 
implies that modulo V{3} the entry is of the form 

d ^ 

fi'\t,0 = i E ^^AtMj + E,=i /^^(*)^'^ + t(^) "i°d r{3} (5.34) 

with aij(t) e T{0}, I3,{t) e T{0} and 7(t) e T{1}. On the other hand, modulo 
0{t^^) the eigenvalues of Fk{t,^) and of A{t,£_) coincide. Since (B3) is a spectral 
condition implying that is a local (quadratic) minimum of an eigenvalue branch 



contained in the complex upper half-plane, some terms in (5.34) have to vanish. In 
particular we see that (3i{t) has to decay, f3i{t) S T{1}, and also that the real part 
of the quadratic matrix (Qfij(i))ij- is positive definite modulo T{1}. The latter is 
a direct consequence of the non-degeneracy of that minimum. 

Corollary 5.3. Assume (Bl), (B2) and B(3) and letk>l. Then modulo 'P{3} the 
upper-left corner entry of F]^{t,S,) satisfies 

/f'^ (t, i) - i^c^m + (iitVi + lit) mod V{?>} (5.35) 

with a{t) e T{0} (8) C"^" having positive definite real part uniform in t > to for 
to sufficiently large, /3(t) G Tjl} «) C" and -/(t) eT{l}. 

Later on we will see that the essential information for deducing diffusion phe- 
nomena and related asymptotic properties is contained in the terms described in 
this corollary. Following methods of [ST] it is possible to relax assumption (B2) 
to the weaker requirement that is simple eigenvalue of B(t) for all t and uni- 
formly separated from the remaining part of the spectrum of B{t). Then a block- 
diagonalisation scheme can be established which separates the corresponding mode 
and proves analogues of the above statements. 

5.3. Asymptotic integration and small frequency expansions 

We consider the transformed problem in V^''^{t,Q — A^fe(t,^)F^°-'(t,^), 

J^tV^''\t,0 ^ {-Dit) + Fk{t,0 + Rk+i{t,0)V^''\t:0, (5.36) 
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and reformulate this as integral equation for its fundamental solution. Denoting 
it by £kit,s,S.), we know that it solves the above equation to matrix initial data 

Let 9fc(i, s, ^) be the fundamental solution to the diagonal system Dt —'D{t) — 
Ffe(i,0-Then 

holds true for t> s uniformly within Zoii(e) for e < c^ sufficiently small. Here, 

Sfc(t,s,0 = exp [ij^ f^\e,Ode\ (5.38) 

gives (for k — 2) the fundamental solution to a parabolic problem and Qk{t,s,S.) 
is exponentially decaying as the fundamental solution of a dissipative system. 
Furthermore, the matrix £k(t,s,£,) satisfies the Volterra integral equation 

ffc(t,s,0 = efc(i,s,o+ / ek{t,0,ORk+iio,O£k{o,s,Od9. (5.39) 



We solve this equation using the Neumann series 

£k{t,s,o^ekit,s,o + J2'' I ^kit,h,ORk+i{ti,0 r ■■■ 

/7_1 J S J S 

ekiti-i,ti,^)Rk+i{ti,^)dU- ■ -dh. (5.40) 

This series converges and can be estimated by 

\\£k{t,s,0\\<expn \\Rk+i{e,0\\de\ . (5.41) 

Based on the estimates for the remainder term i?fc+i(i,^) e 7^{fc + 1}, we even 
obtain uniform convergence within the smaller zone Zon(cfc) n {i|^|'^'"'+^'/^ < S} for 
any constant 5 and \\£k{t, s, (,) — &k{t, s, £_)\\ — > as Ck — > for fixed (5 > as soon 
as we choose fc > 1. 

We can obtain a slightly better estimate based on the uniform invertibility 



of ^fc(i, s, ^) for bounded t|^| in consequence of Corollary 5.3 



Lemma 5.4. Let k > 2 and S > 0. Then the fundamental solution £k{t, s, ^) satisfies 
the uniform bound 

\\£k{t,s,0\\<Ck\^k{t,s,^)\, t>s>to, t|CP<<5, (5.42) 

for some constant Cfe > depending on to, S and k. 
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Proof. Multiplying (5.40) by S^T {t,s,(,) in combination with the uniform bounds 
on 9fc(i, s,f) yields directly 

ii£fe(i,s,oii<|Sfe(t,s,e)i||s^'(i,s,oefc(t,s,e) + ---|| 

< \Ek{t,s,0\exp(c' J^ ||i?fc+i(0,OI|de) , 

whenever t|^P < 5. Furthermore, for k > 2 the remaining integral is uniformly 
bounded on this set. D 

5.4. Lyapunov functionals and parabolic type estimates 

In this section we will partly follow the considerations of Beauchard-Zuazua, [3] , 
and explain how condition (B3) of Section 5.2 allows to derive parabolic type decay 
estimates for solutions to the Cauchy problem (5.161. The construction in [3] was 
inspired by the Lyapunov functionals used by Villani 71 . 



Lemma 5.5. Assume (Bl), (B2), (B3). Then all solutions to (5.161 satisfy the 
point-wise estimate 

\\U{t,S,)r < Ce-^*[«l'||C/o(OlP, K] = |el/(0 ^min{|e|,l}, (5.43) 

in Fourier space with constants C and 7 depending only on the coefficient matrices 
Ak{t) and B{t). 

Sketch of proof. The proof follows essentially [31 Sec. 2.2], we will only explain the 
major steps and necessary modifications to incorporate the time-dependence of 
matrices. As the problem is L^-well-posed and dissipative, it suffices to prove the 
statement only for t > to for sufficiently large <o- 

For a still to be specified selection e — (eg, . . . , Cd-i) of positive reals, €j > 0, 
we consider the Lyapunov functional 



L,[C/;t,C] = ||[/(t,OII 



d-l 



^V-^TTi-t f\ n(i\Aii ^Vl 



min{leU^r'}E^J-I'^(^W^(*'§)'"'^(*'^)'^W^(^'§)'^(*'0)- (5.44) 

The uniform Kalman rank condition (B3) implies that for suitable choices of the 
parameters e and for t > to the two-sided estimate 

J||[/(t,0lP<L,[C/;t,C]<4|lt/(t,e|P (5.45) 

holds true. Therefore, all we have to do is to prove the desired estimate for 
hf[U;t,^] which follows from 

9tL,[t/;i,e]+7[e]'L,[&;t,e] <0 (5.46) 

for suitable 7 and suitably chosen family e. 
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Formally differentiating Le[C/;t,^] with respect to t yields the terms con- 
sidered by [3l Sec. 2.2] together with further terms containing derivatives of the 
coefficient matrices. The latter ones are bounded by 

O{t-^)minm,\^\-'}\\U{t,0r, (5.47) 

which is dominated by 7[^]^L(:[f/;i, ^] whenever i|^| > 1 such that choosing 7 
smaller yields the desired bound. It remains to consider i|^| < 1. Here A(t,£,) 
can be treated as small perturbation of B(t) and the diagonalisation scheme and 



Lemma 5.4 yields the corresponding bound. D 



We will draw a consequence from this statement. It is obtained in combination 
with Holder inequality and the boundcdness properties of Fourier transform. 



Corollary 5.6. Assume (Bl), (B2), (B3). Then all solutions to (5.16) satisfy 

\\U{tr)\\L. <C(l + f)-t(p-i'||C/o||iP (5.48) 

for alll<p<2<q<oo and with r > n{l/p — l/q). 

Proof. For |$| <^ 1 the previous lemma in combination with Sobolev embedding 
theorem yields exponential decay under the imposed regularity. Therefore, it is 
enough to consider bounded S,- Then the estimate is elementary, 

||C/(i,-)||, < \\U{t,-)\\,, < ||e-^*l-l^||.||C/o||p- < C(l + t)-t(^-i)||C/o||p 

provided supp Uo{0 C {C : ICI < 1} and i - ^ - J, - i - i. D 

5.5. A diffusion phenomenon for partially dissipative hyperbolic systems 

Now we will combine the estimates of the previous section with slightly improved 
results obtained from the low-frequency diagonalisation. First we construct a par- 
abolic reference problem, whose fundamental solution is given by 22(i, s,^) and 
afterwards we will explain why and in what sense solutions are asymptotically 
equivalent. 

Following Corollary |5.3| it is reasonable to consider the 'parabolic' probleiqj 

dtw = V ■a{t)Vw + j3{t) ■Vw + i'y{t)w, w{to) ^ wq, (5.49) 

for a scalar- valued unknown function wq. To relate the both problems, we observe 
that the first row of £k{t, s, 0) tends to a hmit as t — ;> 00. This is just a consequence 
of the integrability of Rk+i{t, 0) for k > 2. We use this to define 

Wkis)^ lim ej£k{t,s,0). (5.50) 

It is easy to see that Wk{s) — W2{s) for all k. 



^Note, that ReQ:(i) > in the sense of self-adjoint matrices. We will understand parabolicity in 
this sense. 
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Lemma 5.7. The fundamental solution f/j(t, s,^), k sufficiently large, satisfies the 
estimate 

\\£k{t,s,0-^k{t,s,0eiW2{s)\\<Ck{l + t)--2, t>s>to, (5.51) 

uniformly on \^\ < 1. 

Proof. We make use of a constant S > 0, to be fixed later on, to decompose the 
extended phase space into several zones. 

1. If t|^p > i51og< with S chosen large enough, both terms can be estimated 
separately by exp{—jt\^\'^) for some constant 7. This follows for the first one by 



Lemma 5.5 



and for the second by the parabolicity of (5.49) in consequence of 



Corollary |5.3| But then 

e-7t|?P < e-ys log t ^ ^-^5 <t-i^ 7<5 > ^ . (5.52) 

2. If t|^p < 6 for some S, we use the results from the asymptotic integration of 
the diagonalised system. First, we claim that 'E,'^^{t,tQ,£_)eJSk{t,tQ,^) converges 
locally uniform in 1^ as t — ?► ij for t^\£,\'^ — 5. To see this, we use the Neumann series 
representation (5.40) multiplied by the uniformly bounded S^r^(i, s,C) combined 



with Cauchy criterion. We denote the resulting limit as Wfc(s,^) and observe that 
it coincides with Wk{s) for ^ = 0, 

\\Wk{s,0-Wk{s)\\<n, (5.53) 

and satisfies 

\\ejEk{t,s,S,) -^k{t,s,OWk{s,m <t'^- (5.54) 

The first of these estimates follows as uniform limit for estimates of the difference 
S^^(i, 5,0)6]'^ £'fc(i,s,0) — S^^(i,s,^)e7ffc(t, s,^). Indeed, using the Neumann series 
we see that the first terms are equal, the second terms are reduced to the estimate 
\\Rk+i{t, 0) - Rk+i{t, Oil < \£.\{t^^ + ICI*") following directly from the definition of 
the V{k + l}-classes and, therefore, 

||S-i(t,s,0)e7f,(i,5,0)-S-i(t,s,C)e7ffc(i,s,OII 

<iei l\e-'' + \^\'')de+... (5.55) 



and the right-hand side is uniformly bounded by |^|. Taking limits proves the 
estimate. The second estimate is similar. Again using the Neumann series we see 
that this difference can be estimated by 

\\^l\t,s,Oel£k{t,s,0-Wk{s,m 



5 



< r \\Rk+i{T,0\\expU \\Rk+i{0,m<iO]dT<t-i (5.56) 
due to ||i?fe+i(t,e)|| ^ 0{t-^ for fc > 2 and t\C\^ < S. 
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Combining both of the above estimates and using that the other rows in £k 
are exponentially decaying we get 

<i-5+e-^*l«l'|^| <t-3. (5.57) 

3. It remains to consider the logarithmic gap between both parts, i.e., 5 < t\^\'^ < 
i51ogi. Here we use that for k sufficiently large the remainder term Rk+i{t,£,) 
decays as t~'^~^+'^, while the polynomial growth rate of S^^(i,^) is independent 
of k for large k. Choosing k large enough, the Neumann series argument gives 

\\E-\t,s,OeJSk{t,s,o - Wk{s,m < I ' m-\9,3,ORk+i{o,0¥e < r'^ 

(5.58) 

with t^ defined by i^\£,\'^ = (Jlogi^ and Wk{s,(,) ~ lim^^j.^ ^'f^^{t,s,(,)elEk{t,s,S,)- 

The existence of the latter limit follows for large k and again Wk{s,£,) — Wk{s) 
coincides up to order 0(1^1). D 

To obtain a statement in terms of the original equation, we introduce K{t, S,) — 
M{t)N2{t,£)ei. By definition we have K{t,£,) £ 'P{2}. We define further wq = 
WUq in such a way that we cancel the main term of the solution within Zcii{ck) H 
{t\^\ < S}, i.e., we define 

wo = W2ito)N^\to,OM-'ito)£{to,0,OxiOUo (5.59) 

with xiO e Cf^iW), xiO = 1 near f = and suppx C ^^^(O). Then the estimate 
of Lemma |5.7| implies the following statement. The logarithmic term is caused by 
comparing Ek(t,s,£,) with S2(i, s,^). 



Corollary 5.8. Let U{t,x) be solution to (5.161. The the solution w{t,x) to (5.49) 
with data given by ( |5.59[ ) satisfies 

\\U{t, •) - K{t, D)w{t, ■)\\2< C'{1 + i)-3 log(e + 1). (5.60) 

As only the upper-left corner of the diagonalised problem is of interest, As- 
sumption (B2) can be relaxed and it is sufhcient to guarantee block-diagonalis- 
ability. We will not go into detail here, but remark that it is sufficient to have that 
the eigenvalue is uniformly separated from the remaining spectrum of the family 
B{t). Then the diagonalisation scheme can be modified based on the method of 
|67j and an analogous result obtained. 
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6. Examples and counter-examples 

Both, in Section H and Section [5] we made symbol like assumptions on coefficients, 
e.g., we considered hyperbolic systems 

n 

D,C/ = ^A,(i)D,,t/ (6.1) 

fc=i 

for coefficient matrices Ak{t) G T{0}, meaning that derivatives of the coefficients 
are controlled by 

mA,it)\\<C,(^^^ . (6.2) 

We will use this section to show that assumptions controlling the amount of oscilla- 
tions in the time-behaviour are in fact necessary in order to control the large-time 
behaviour of the energy. Yagdjian pointed out in [70] the deteriorating effect time- 
periodic propagation speeds might have on energy estimates and more generally 
on the global existence of small data solutions for nonlinear wave models. The key 
idea behind these is based on Floquet's theory for periodic ordinary differential 
equations and in particular Borg's theorem. 

We will recall these first in connection with a very simple model and then 
give an idea how these can be used to prove sharpness of estimates and sharpness 
of assumptions on coefficients for wave models. 

6.1. Parametric resonance phenomena 

We will restrict ourselves to the simple model 

utt — a {t)Au = 0, ii(0, •) = lio, ut{0, ■) — ui, (6.3) 

of a wave equation with variable propagation speed. Basic assumption will be that 
a^{t) is positive, smooth, periodic 

a{t + l)=a{t), (6.4) 

and non-constant. Using a partial Fourier transform with respect to the 2;-variable, 
this is equivalent to Hill's equation 

utt + \^\^a\t)u^0 (6.5) 

with spectral parameter |^p or, equivalently, the first order system 

for U — {\i\u, Djm)^. Its fundamental matrix is again denoted as £{t, s, ^). Due to 
periodicity, it is of interest to consider the monodromy matrix M{£,) — £{1,0,£,). 



Since the trace of the matrix in (6.6) is zero and hence det7V{(^) — 1, it has 



either two non-zero eigenvalues of the form exp(±K(^)) or is a Jordan matrix to 



the eigenvalue 1. The number k(^) is called the Floquet exponent of (6.5). Its 



importance stems from the following lemma. It can be found in a similar form in 
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Lemma 6.1 (Floquet-Lemma). Assume that M{£,) is diagonalisahle with Floquet 



exponent k{^). Then (6.5) has a fundamental system of solutions of the form 

e±*^(«)/±(t,C) (6.7) 

for non-vanishing 1-periodic functions f±{t,^) depending analytically on ^. 



Proof. We will give a proof based on the system (6.6). The same method applies 
to arbitrary periodic systems of differential equations. By assumption we find an 
invertible matrix N{^) such that 

M{(,)N{0 = N{£_) diag(e'"(«\ -e''^)) = N{^) exp {K{(,)) (6.8) 

holds true with K{£^) = diag(K(^), — k(^)). We use the fundamental matrix to 
define the function 

F{t, = £{t, 0, 0N{0 exp {-tKiO) N~\^). (6.9) 

Then a simple calculation shows that the periodicity implies £'(i + f,l,^) = 
£{t,0,£,) and, therefore, we obtain 

F{t + l,0 = £(t,0,aM(OA^(C)e-^(«^e-*^(«)iV-i(O = F{t,0 (6.10) 



based on (6.8 ). Furthermore, by construction we see that any solution of the system 
(6.6) is of the form 

U{t, = Sit, 0, OUo = F{t, exp i-tKiO) Uq. (6.11) 

The statement of Lemma |6 . 1 1 follows by looking at individual entries of this matrix 
in combination with the definition of U. D 

Since a{t) was assumed to be real, the above representation implies a sym- 
metry for the Floquet exponents. Either, they are both imaginary or they are 



both real. If they are imaginary, all solutions to (6.5) remain bounded, while for 



real Floquet exponents an exponential dichotomy appears. Looking at the original 



equation (6.3) this means that for appropriately chosen initial data, solutions will 
have an exponentially increasing energy. 

For us the following statement is of interest. In the present form it is due to 
Colombini-Spagnolo fB^, for more detailed results on Hill's equation see Magnus- 
Winkler [27,. 

Lemma 6.2 (Borg's theorem). Assume the coefficient a{t) is locally integrable, 1- 
periodic and non- constant. Then there exists an open interval J C M+, such that 
for all 1^1 G J^ the Floquet exponent satisfies k(^) > 0. 

In combination, the above lemmata imply that we can find initial data such 
that the solution is exponentially increasing. Note, that no assumption on the size 
of the coefficient was made, we only used that a{t) is not constant. 

Corollary 6.3. There exist inital data uq, ui G ^(M") such that the solution u{t, x) 



to (6.3) satisfies 



liminf ^ ^ ' ^ > 0. (6.12) 
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Figure 3. Sketch of the constructed coefficient a{t). 



6.2. Construction of coefficients and initial data 



We stay with the model (6.3), 

utt - a'^it)Au = 0, u{0,-)^uo, Mt(0,-)=ui, (6.13) 

but drop the periodicity assumption. The above sketched instability mechanism 
will be used to construct coefficient functions and data for prescribed energetic 
behaviour. 

To be precise, let r^ /* cx), Sk, »7fe and Uk be sequences such that 

Tk+i > Tk + Sk, Vk < 1, nk £ N>o. (6.14) 

The sequences t^ and 6k define disjoint intervals Ij, — [Tk,Tk + 5k]. We will let our 
coefficient oscillate at least Uk times on the interval Ik ■ Let therefore <j) G C^ (M) 
be supported in the interval (0, 1) and non- vanishing with |</>(i)| < 1 and denote 
by b{t) its periodisation, b{t) = (j){t — [tj). Then we consider 



«(*) = 
The function a{t) is bounded, 



Vkb{'jf{t-rk)), 



telk, k 

otherwise. 



1,2,. 



< a{t) < 2, 



and satisfies the estimate 



^nk-^l 



Ok' 



teik 



(6.15) 

(6.16) 
(6.17) 



the constants C^ are independent of k. 

The freedom in the choice of sequences can be used to construct examples of 
coefficient functions in order to prove sharpness of energy estimates. For instance, 
we obtain a{t) e T{0} for the choice Tk '^ k, Sk ''^ k/2, rjk = 1, and n^ = 1 while 
increasing Uk to nk = [/c^^"], a £ (0, 1), yields the weaker estimate 



Ha{t)\ < C, 



1 



(6.18) 
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If we want to show lower bounds on the energy or the sharpness of energy 
estimates, we have to construct initial data. One particularly simple idea is to use a 
sequence of initial data Mo,fcj uu^ £ J?*(M") which leads to an exponential increase 
of the energy within the interval Tk and still satisfies good estimates in the earlier 
intervals Ik', k' < k. The latter can be achieved by adjusting the sequences in the 
definition of the coefficient function. This idea was employed in [T2], [Ej to discuss 
the sharpness of the results. 

The following statement is entirely in terms of the defining sequences and 



provides such a lower bound on the energy behaviour of solutions to (6.13). We 
use r/fc = 1 for simplicity in order to reduce the statement to Borg's theorem in 
the form of Lemma 16.21 



Lemma 6.4. Assume a{t) is defined in terms of the above sequences r^, 6k and 
Uk- Then there exists a sequences of initial data uo.fcj ui.k with normalised initial 



energy E{uk',0) — 1 such that the solution to (6.13) satisfies 

fc-i 
logE{uk,Tk + Sk)>2Knk-2c'^ni (6.19) 

i=i 

with K > small and c — sup^ ' AV' > 0. 

Proof. We first outline the main strategy. We choose the sequence Uk in such a way 
that Uk{Tk,^) restricted to the interval Xk is an exponentially increasing (eigen-) 
solution to the periodic problem on that interval. To relate different fc, we note 
that all appearing periodic problems are rescaled versions of each other and we use 
the same instability interval J^ for all of them. We follow these solutions backward 
to the initial line and adjust them in such a way that they are normalised in initial 
energy. 

To be precise, let J^ be an interval of instability of the periodic problem 

vtt + \^\\l + bit)fd ^ 0. (6.20) 



Let further v{t,^) be an exponentially increasing eigensolution of the form (6.7) 
supported inside the interval J^ such that on supp v the Floquet exponent satisfies 
/^(O > K > for some positive constant k. We can rescale v to 

Uk{t,0=l^kvC^{t~Tk),—0, ?7fcGK+, (6.21) 

Ok rik 



such that Uk{t, x) solves (6.13 1 on Ik- We will use ^k to normalise solutions. Equa- 
tion ( 6.21[ ) implies 

nuk,Tk + 5k) = c2«-"E(ufe,Tfe). (6.22) 

Outside of IJ^If the energy is conserved, while it might decrease or increase on 
the intervals Ti. Therefore we obtain the lower bound by looking at the worst type 
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behaviour as estimated by Gronwall inequality from 

\a'it)\ 



dt / (|wt|' + a'W|Vw|2)dx 



= 2- 



a{t) 



a\t)\S/uY<lx 



<2^# / {\ut\^ + a\t)\^u\^)Ax, (6.23) 



a{t) 



W(t)\ _ «ic„T. I'''(^)l 



i.e., E(ufc,T£ + 5t)>e '^^"''E{uk,T() with q = sup^gx, ^^(tf = "57 sup^ -^ 
combination this yields 



. In 



fe-i 



logE(Mfe,rfe) -logE(ufe,0) > -2c^ni 

1=1 



and the statement follows by combining (6.24) with (6.22 1 



(6.24) 

D 



When using this construction one has to pay attention to some particular 
facts. First, the bad behaviour of solutions is localised to |^| ~ n^jS-k — >■ on 
T-k = \rk^Tk + (5fe] — >■ c». In view of (6.17) this is related to the decomposition of 



the phase space into zones and the bad increase in energy happens close to the 
boundary of the hyperbolic zone. Second, the choice of sequences does matter. 
Even if it can not be seen in the statement of the previous lemma, one wants to 
construct coefficient functions violating conditions as closely as possible in order 
to draw interesting conclusions. 



We draw one consequence from the above lemma. When considering (6.13) 



with a(t) e T{0}, then solutions satisfy the a global generalised energy conserva- 
tion property in the sense that 

E(u;i) «E(m;0) (6.25) 

uniform in t and with constants depending only on the coefhcient function a(t). 



This follows from Theorem 4.6 and was originally proven in Reissig-Smith [46]. If 
we use the sequences rj, = ct*^, bk = o^~^ and n^ = \u'^^'\ for some g > and a 
given parameter ct > 0, we obtain a coefhcient a(t) satisfying the estimate 



|Dfa(i)| 



teifc 



<C, 



\ai^^ 



Cp 



l + t 



[l~q)t 



The solution to this problem can not satisfy (6.25), as 



qk 



fc-i 



^na'i^ 



' cr9- 1 



oo, 



(6.26) 



(6.27) 



provided a is chosen large enough in comparison to k and c. 

A slight modification of the argument allows to show that the assumption 



a{t) e 7I/{0} for some u > Q is also not sufficient to deduce the estimate (6.25) 
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7. Related topics 

Most of the result presented here were based on diagonahsation procedures in 
order to deduce asymptotic information on the representations of solutions. This 
is natural and has a long history in the study of hyperbolic equations and coupled 
systems. For diagonahsation schemes in broader sense and their application we 
also refer to [3T]. Some more applications are discussed there too. 

Our main concern was the derivation of energy and dispersive type estimates 
describing the asymptotic behaviour of solutions to hyperbolic equations. As it 
is impossible to cover all directions appearing there, we refer also to the expos- 
itory article [68] for a discussion on current results about energy type estimates 
for wave models with bounded coefhcients. Wave models with unbounded coeffi- 
cients (meaning polynomially increasing propagation speed or even exponentially 
increasing propagation speed) have been extensively studied by Reissig [40] and 
Reissig-Yagdjian, 01], I^, [S], [H], see also Galstian [TD]. 

There is an interesting duality to be observed here. Wave models with in- 
creasing coefhcients and their large-time behaviour are intimately connected to 
the well-posedness issues of weakly hyperbolic equations. In a similar way, models 
with bounded coefficients are related to well-posedness statements of equations 
with log-Lipschitz coefhcients and their microlocal analysis. See, e.g., the results 
of Kubo-Reissig [24^ or Kinoshita-Reissig ^5j where the key ingredient of the 
consideration is a diagonahsation scheme for pseudo-differential versions of the 
symbol classes introduced in Section [44| applied locally in time. 

The considerations of the generalised energy conservation property in Sec- 
were inspired by the work of Hirosawa, '12\ and later joint work with 



4.5 



tion 

him. |13]. |14j . Main question arising there is to what extent symbolic conditions 
on coefficients can be weakened without loosing uniform bounds on the energy. 
As the presented results are sharp, weaker estimates for derivatives have to be 
compensated by additional stabilisation conditions. There has been recent work in 
this direction by D'Abbicco-Reissig [2 for 2 x 2 hyperbolic systems generalising 
results of 14 and also ongoing joint research of the second author with Hirosawa. 
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